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ABSTRACT 


The  guiding  of  waves  in  interferometrically  produced  waveguides  is  in- 
vestigated. A very  general  technique  of  analysis  was  developed.  Computer 
calculations  based  on  this  approach  were  used  to  study  leakage  from  guiding 
channels  and  the  effect  of  pulse  shape.  Experimental  devices  were  construc- 
ted. A new  concept,  modular  optics,  was  conceived. 


SUMMARY  AND  ACCOMPLISHMENTS 

During  the  course  of  this  contract,  we  pursued  the  following  objectives, 
as  stated  in  our  plan  of  work  in  our  initial  proposal: 

1.  We  will  adapt  the  techniques  of  volume  holography  to  the  devel- 
opment of  multichannel  optical  waveguiding  structures.  We  will 
expieriment  with  a number  of  the  most  successful  materials,  such 
as  dichromated  gelatin,  thick  plastics  (e.g.,  polymethyl  metha- 
crylate, and  various  photopolymers, 

2.  We  will  develop  methods  for  controlling  the  characteristics  of 
these  guiding  channels  through  appropriate  modulations  of  the 
interfacing  beams.  In  particular,  we  will  form  cylindrical 
lenses  within  these  channels,  and  will  explore  methods  to  per- 
mit these  lenses  to  have  focal  lengths  different  for  each  chan- 
nel. 

5.  We  will  develop  methods  for  efficient  coupling  of  light  beams 
into  these  channels. 

. On  the  basis  of  what  our  studies  indicate  is  realiziable,  we 
will  explore  the  feasibility  of  various  forms  of  optical  pro- 
cessing devices. 

Gome  of  these  proved  quite  attainable,  whereas  others  were  found  to  be 
difficult  and  success  was  only  partial.  In  addition,  various  new  ideas,  net 
anticipated  in  the  original  statement,  but  broadly  falling  within  its  scope, 
were  conceived. 

We  list  our  various  accomplishments,  as  they  relate  to  the  four  items  of 
our  work  plan. 

Our  principle  efl’ort  was  directed  toward  item  1.  Our  investigation  of 
applying  the  techniques  of  volume  holography  to  the  development  of  multichannel 
guiding  structures  was  done  mostly  on  an  analysis  basis. 

A major  problem,  recognized  at  the  outset,  was  to  find  a theory  to  aid  in 
the  annlysir.  of  guiding  in  closely  spaced  channels  having  index  profiles  of 
the  sort  attainable  with  ini.erferometric  construction.  Gince  the  channels  are 
to  be  constructed  by  recording  interference  fringes,  and  since  the  intensity 
j)rofile  is  nor-malLy  sinusoidal,  we  expect  that  the  resulting  channels  will  ex- 
hibit a sinusoidal  index  profile.  The  analysis  of  the  waveguiding  properties 
of  such  channels  is  not  easy,  and  we  felt  that  the  conventional  means  were  not 
very  satisfactory.  To  this  end,  we  developeii  a new  anii  essentially  differer.t 
approach,  which  is  highly  versatile  and  quite  simple  to  apply.  This  is  the 
thin  grating  decomposition  technique  (hereafter,  TGD)  developed  originally  for 
analysis  of  li  i f frac+.i  ■ n from  volume  gratings,  but  readily  adapted  to  the  multi- 
channel waveguide  problem.  This  technique,  as  it  applies  to  waveguiding.  is 
described  in  appended  paper  No.  1,  by  Leith,  ^ ^1. 

We  desire(i  a theory  which  would,  without  more  than  a change  of  parameters. 


L^di 
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describe  both  Bragg  diffraction  from  thin  gratings,  and  the  guiding  effects 
that  occur  with  thick  materials  and  high  index  modulations.  We  desired  a 
theory  that  would  describe  the  transition  region,  when  partial  guiding  ef- 
fects are  observed,  and  which  would  also  describe  all  the  effects  of  the 
various  elements  that  we  wish  to  introduce  into  the  guiding  channels. 

The  TGD  method  fulfills  all  of  these  requirements,  and  we  have  used  it 
extensively . 

In  our  first  report  (Jieport  No.  1,  Annual  Report,  April  197^  - March  19'^5) 
we  used  the  TGD  method  to  examine  leakage  between  channels  and  pulse  shaping. 
The  bulk  of  this  work  was  published  in  Applied  Optics ; the  paper  is  appended 
as  paper  No.  2. 

This  work  has  been  continued,  resulting  in  another  paper,  submitted  for 
publication  in  Applied  Optics  (appended  paper  No.  5).  We  summarize  the  re- 
sults of  that  paper  as  follows:  We  assumed  a sinusoidal  index  profile 

throughout,  and  we  generally  assumed  the  input  energy'  to  have  a shape,  in  the 
direction  transverse  to  the  channel,  of  one  cycle  of  a sinusoid.  We  found 
that  ■t.he  index  modulation  required  of  the  guiding  structure  to  rise  rapidly 
with  decreasing  channel  spacing,  being  roughly  quadratic  in  the  region  of  90- 
190  channel/mm,  and  to  rise  m.uch  more  rapidly  above  150/mm.  so  that  as  a 
practical  matter  channel  densities  above  about  200/mm  are  not  feasible. 

We  examined  the  guiding  characteristics  of  various  pulse  profiles,  cthe;’ 
than  the  sinusoidal  one  previously  noted.  One  interesting  observation  is 
that,  regardless  of  the  spatial  width  or  shape  of  the  input  pulse,  the  out- 
put pulse  tended  to  converge  to  one  of  sinusoidal  profile  matching  in  scale  to 
that  ol’  t.he  guiding  channel,  that  is,  if  the  channel  profile  is  of  the  form 
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t(v)  = COS  dn  ^ 

This  i nt,er’f-st i ng  phenomenon  can  be  explained  in  terms  of  the  basic  mode  equa- 
tions for  the  sinusoidally  profiled  waveguide. 

i'hese  are  a few  examt  les  of  results  from  our  extensive  study  of  guiding 
cha  racter  i st  i cs  using  the  T f metho'i. 

Tn  aiidition  to  the  theoretical  study,  we  carried  out  a continuous  experi- 
mental investigation  of  materials.  We  fi rst  used  polymethyl  methacrylate  an  1 
various  photopolymers , tlien  tried  li ichromai.ed  gelatin.  The  latter  is  (iesir- 
able  because  of  the  large  reCractive  index  changes  that  can  be  induced  in  it. 


at  least  10  X that  inducible  in  polymethyl  methacrylate.  Dichromated  s^elatln, 
hi>wever,  has  its  own  special  problem,  the  principle  one  being  that  samples  of 
thickness  greater  than  ‘)0  to  100  microns  are  not  feasible. 

We  first  explored  commercial  plates  which  came  in  thicknesses  of  15  to  JO 
microns.  These  were  Kodak  6U9  photographic  plates,  and  we  prepared  them  for 
our  purpose  by  removing  the  silver  and  sensitizing  the  gelatin  with  ammonium 
dichi'omate.  This  is  a standard  procedure.  We  then  proceeded  to  prepare  our 
own  plates  by  mixing  our  own  emulsion  and  applying  it  to  glass  plates.  We 
could  thus  make  gelatin  films  of  any  thickness  we  chose.  This  became  a com.- 
plicated  experimental  program;  we  had  to  cope  with  many  problems,  such  as  bend- 
ing of  the  gelatin  to  the  substrate;  drying  of  the  gelatin  without  cracking  or 
peeling,  sensitizing  the  gelatin  uniformly  throughout  its  depth,  and  prepara- 
tion of  specimens  which  would  remain  bonded  to  the  glass  substrate  throughout 
the  entire  procedure,  including  sensitization,  development,  washing,  and  dry- 
ing. 

Our  goals  were  to  try  to  form  dichromated  gelatin  plates  as  thick  as  pos- 
sible, while  still  achieving  index  modulation  sufficient  to  produce  waveguiiies. 
We  succeeded  in  achieving  thicknesses  of  about  50  microns  while  still  achieving 
the  required  index  modulation.  This  thickness  is  adequate  for  some  of  cur  pur- 
poses, but  not  for  others. 

The  advancement  of  dichromatefi  gelatin  technology  is  a long-term  project 
which  is  important  to  many  of  the  activities  of  our  laboratory.  We  expect  cur 
research  efforts  in  this  area  to  continue,  with  progress  being  rather  gradual. 
Similarly,  progress  elsewhere  is  continuing,  with  gradual  improvement  in  the 
technology . 

Our  effort  on  item  2 of  the  work  plan  was  also  extensive,  although  re- 
stricted to  single  channels. 

Principally,  we  developed  the  concept  of  modular  optics,  which,  while 
falling  within  the  scc'pe  of  that  item,  is  we  believe  a worthwhile  concept  in 
its  own  right  and  goes  beyond  the  scope  of  the  work  statement.  This  concept, 
stated  briefly,  is  that  one  may  form  (in  our  case,  interferometrically ) opti- 
cal element.s  of  various  sorts  which  can  be  attached  by  index-matching  fluid 
coupling  to  a slab  wavegui'ie  whereupon  a guidel  wave  is  pulleri  out  of  the  slab, 
is  operated  upon  by  the  modular  element,  and  then  in  its  modified  form  is 
coupleii  back  into  the  waveguide. 

The  mo  iular  element  can  be  a bcim  splitter,  a lens,  a mode  converter,  or 
any  other  element  that  can  be  fabricated  as  an  integrated  optics  structure. 

In  this  manner  optical  systems  can  be  constructed  that  are  in  fact  integrated 
optic.s  systems,  having  the  ruggedness,  stability,  and  compactness  of  such  sys- 
tems, but  also  having  the  adjustability  of  conventional  optical  systems. 

This  concept  has  beer,  disclosed  in  a patent  applicati  'n,  and  at  present, 
four  papers  have  been  published  or  submitter!  on  various  aspects  of  this  con- 
ceid..  To  ilate,  we  have  constructed  a m 'tular  optics  interferc'meter , and  a 
double  mo<ie  coupler. 


On  the  third  item  of  the  work  plan,  the  coupling  into  the  guides,  we  made 
considerable  progress.  Of  course,  the  modular  optics  devices  just  cited,  come 
under  this  item  also.  In  addition,  we  conceived  a technique  for  producing  a 
high  efficiency  grating  coupler  within  the  waveguide.  We  analyzed  this  struc- 
ture using  the  TGD  method,  and  found  from  computer  analyses  that  the  coupling 
efficiency  could  be  high,  nearly  100^.  This  work  is  described  in  the  appended 
paper  by  Wang  & DiLaura,  submitted  to  Applied  Optics . 

At  the  basis  of  all  of  our  analysis  is  the  TGD  methofi . Developed  by  Alfer- 
ness  to  analyze  third  transmission  gratings,  the  scope  of  its  applicability  has 
subsequently  been  enlarged,  having  been  used  to  analyze  grating  couplers,  re- 
flection gratings,  diffraction  at  the  third  Bragg  angle  of  non-linear  grating, 
and  distributed  feedback  lasers  where  the  modulation  index  of  the  active  me- 
dium varies  sinusoidally  rather  than  in  discrete  steps. 

Finally,  much  of  our  device  analysis  falls  within  the  scope  f item.  1. 
the  final  item,  of  our  work  plan. 

A total  of  11  papers  were  proposed  for  publication  from  work  supporte  ; 
either  entirely  or  in  part  by  this  grant.  Five  of  these  are  now  in  print  and 
the  remaining  six  are  pending.  In  addition,  two  doctoral  ii ssertations  (G. 

Kung  and  W.  Wang)  were  supported  by  this  work.  The  papers  are  presenteri  in 
the  Appendix. 

The  principle  obstacle  to  the  fullest  realization  of  our  objectives  is 
the  materials  problem;  materials  having  the  ideal  characteristics  (i.e.,  very 
thick,  high  transmittance  before  development  for  the  exposed  stage)  and  high 
index  modulation  are  not  available.  We  therefore  (a)  worked  within  the  mate- 
rials limitation;  (b)  performed,  on  a modest  scale,  materials  development;  and 
(c)  carried  out  analysis  and  design  on  the  basis  of  material  characteristics 
not  yet  available.  In  addition,  some  of  our  results  (i.e.,  the  TGD  analysis 
and  the  modular  optics  technique)  have,  we  believe,  import  well  beyond  the 
scope  of  our  grant  objectives. 
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The  ihin-praliny  decoinpusiiion  melhod  ol  thick  ttraling  analysis  is  used  to  analyze  puidiny  elTects  in  thick  phase 
piatinps  having  sinusoidal  retractive  index  modulation. 


Intert'eromctrically-produced  phase  gratings  tan 
be  used  as  wave-guiding  structures,  as  demonstrated 
by  Rosenberg  and  Chandross  1 1 1.  Rosenberg  gives  an 
analysis  of  this  phenomenon  by  ray  tracing  methods 
[2|.  We  note  tliat  the  structure  that  produces  Bragg 
diffraction  is  the  same  one  that  produces  guiding,  ex- 
cept that  for  guiding  to  be  manifested,  both  the  thick- 
ness and  index  modulation  must  be  relatively  large. 

We  also  note  that  the  basic  Bragg  diffraction  equation 
and  the  basic  equation  describing  guiding  in  a slab 
waveguide  are  similar.  We  have  therefore  sought  a 
single  analysis,  preferably  a grating  analysis  based  on 
physical  optics,  that  describes  both  phenomena,  and 
that  will  give  insight  into  the  transition  region,  when, 
with  increasing  refractive  index  modulation,  the  Bragg 
diffraction  process  produces  con nnement. 

In  the  diffraction  process,  a plane  wave  impinges 
on  the  thick  grating  at  the  Bragg  angle  (fig.  1 ).  When 
the  wave  enters  the  structure,  a diffracted  wave  is 
formed;  as  the  two  waves  propagate  througli  the 
structure,  there  is  a continuing  exchange  of  energy 
between  them.  Continuing  with  this  heuristic  ap- 
proach, we  expect  that  when  most  or  all  of  the  ener- 
gy in  a plane  wave  is  rellected  from  the  first  Bragg 
plane  encountered,  energy  initially  contained  between 
two  Bragg  planes  is  confined  between  these  planes 
until  the  wave  emerges  from  the  structure.  Unfortu- 
nately, the  usual  analyses  for  dilfraction  from  thick 
gratings  assume  uniform  plane  waves.  The  analyses 
predict  the  intensity  ol  the  emerging  waves,  but  give 
no  insiglii  into  the  paths  that  the  ligjtl  traveled.  In 
addition,  the  closed-form  solutions  are  restrictive; 


I ip.  I.  Sinusoidal  thick  pratinp,  showinp  incidc’nt  and  dil- 
Iraclcd  waves. 
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I ip.  2.  Spatial  lipht  pulse  incident  on  pratinp  surface. 

Kogelnik's  coupled  wave  analysis  |3).  for  example, 
allows  for  only  two  plane  waves,  and  these  must  each 
be  at  or  near  the  Bragg  angles.  We  require  an  analysis 
that  enables  us  to  propagate  a small  spot  of  light  (a 
spatial  pulse)  througli  the  structure. 

Alferness  (4|  has  described  an  extremely  general 
method,  which  allows  the  plane  waves  to  be  incident 
at  any  angle.  It  is,  however,  a plane-wave  analysis, 
but  it  can  be  modified  for  oui  requirement.  This 
method  is  called  the  thin  grating  decomposition  meth- 
od, and  is  conceptually  quite  simple,  yet  in  all  tests 
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iIkiI  Ikivo  Iil'cii  iii.Kk’.  Iiulli  liy  wuy  ul  c\|>oi iiiii.'iil 
and  by  compaiisDii  with  S|X’cial  cases  llial  can  be 
ajialysed  by  oilier  nielhods,  the  Allerness  nielliod 
has  proved  accurale.  The  inelliod  involves  deconi- 
|soslng  a ihick  graliii)’  inlo  a slack  ol  thin  gratings, 
each  id  which  is  siilTiciently  tliin  that  thin  grating 
theory  holds,  l ach  tliin  grating  dilTracts  the  incident 
plane  waves  inlo  various  ordcis,  depending  on  the 
grating  sirucliire.  These  dilTracted  waves  are  the  in- 
put lor  the  next  thin  layer.  Typically,  the  decom- 
position involves  about  30  to  100  thin  gi  atings;  thus, 
the  analysis  rei|uires  a computer. 

In  our  procedure  (tig.  2),  a spatial  pulse  impinges 
on  the  source-side  surface  of  the  grating.  The  pulse 
is  l ourier-deconiposed  inlo  a spectrum  ol  plane 
waves  at  various  angles.  The  Allerness  theory  per- 
mits us  to  propagate  each  one  through  the  strucluie. 
and  to  combine  all  ol  the  emerging  waves,  both 
direct  and  dilliacled,  to  Ihid  the  lesullanl  intensity 
disli  ibulion. 

Since  lor  calculalional  purposes  we  ate  restiicled 
to  a Iniile  nunibei  of  plane  waves,  we  should  choose 
an  aniplilude  piollle  lor  the  pulse  which  is  not  only 
realistic,  but  which  can  also  be  described  without 
an  excessive  nuniber  ol  components.  A suitable 
wavelorni  is 

,/ (.V ) - ( i j cos  2 nx  J)  reel  .v/</  , ( 1 ) 

where  J is  the  grating  period,  and  reel  x/J  = I for 
|,v|  < \ J.  and  IS  zero  otherwise.  This  waveform  has 
the  advantage  that  not  only  is  it  conrined  between 
the  limits  |.v|  = ^ J.  but  also  that  its  spatial-frequen- 
cy spectrum 

) = \ J sine  JJ\.  + 4 (/  sine  + I UJ) 

+ 4 </  sine  </(./^  l i/  I (2! 

IS  very  nearly  conlined  between  the  values  = 2IJ. 
(Note  that  sine  v (sin  rr.vl/rr.rl. 

We  sample  this  s(X’clruni  at  intervals  I TV.  giving 
4/’  I samples  between  the  limits  l 2/T  (tig.  3). 
giving  a sampled  s|K’clrnni 

) = f'Uy  > comb  /V/^  . (3) 

wheie  comb  aj^  = li  6(_/^  ii  i/l  and  ft  is  the  Diiac 
della  function.  The  sampling  lesulls  in  a niodilled 
input  Innclion 
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I i^.  3.  Spulial  Ircquencj  spA*tlrum  of  pulse.  Only  the 
envelope  is  shown.  63  discrete  trequcncies  reside  under  the 
envelope. 


/,,(.v)  = |(j  +5  cos  2-nx/J)  lect  x/<J\  * comb  x/PJ 
2P 

= S a, , cos  2niix/J  . (4) 

2P 

This  function,  being  composed  of  a finite  nuniber  ol 
harmonically-related  spatial  frequencies,  thus  physi- 
cally corresponds  to  a set  ol  4/^  1 plane  waves,  with 

incidence  angles  0 = sin  ' iiX/V.  These  we  propaga'e 
througli  the  grating  structure,  using  the  .Allerness 
method.  Some  of  these  components  satisfy  nearly  oi 
exactly  the  Bragg  condition  and  produce  dilfracted 
waves  at  the  output ; others  do  not.  At  the  output 
plane,  or  exiting  surface,  .ill  incident  componeiils. 
along  with  the  new  components  generated  by  difli ac- 
tion. are  combined,  thus  forming  the  output  wave- 
form. 

Some  precautions  should  be  observed  in  designing 
the  experiment.  The  sampling  process  forms  a periodic 
function  from  the  single  pulse  J : thus,  the  input  is 
a sequence  of  pulses  of  the  form  ,/,  separated  by  inter- 
vals PJ.  We  should  choose  /’sulficiently  large,  i.e  . 
take  sufficient  samples,  that,  over  the  propagation 
distance,  the  pulse  spreading  due  to  diffraction  does 
not  result  in  an  overlap  ot  adjacent  pulses.  The  giealei 
the  nuniber  ol  samples.  4/’  I.  the  greater  the  separa- 
tion distance  between  pulses. 

The  calculated  results  aie  shown  in  fig.  4.  We  have 
assumed  a giating  with  a sinusoidal  variation  of  re- 
fiactive  index,  with  index  variation  along  the.v  diiec- 
lion  only  (the  unslanied  liinge  case).  The  |X’rtinenl 
parameters  are:  giating  spacing  T - I '100  mm.  grating 
thickness  = 0.54  mm.  wavelength  X = 0328  ,\.  and  a 
refractive  index  of  1.50.  The  progiam  broke  this  into 
18  thin  giatings  of  thickness  0.03  mm;  such  gratings 
satisfy  Kogelnik's  y-factoi  ciiteiion  for  a thin  giating 
(3|.  The  output  pulse  is  shown  in  a to  d tor  various 
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I'ig.  4.  Intensity  distribution  at  output  surface,  for  various  values  of  A/r:  (a)  An  - 0 (no  pratinp):  (b)  An  = 10"^:  (c)  An  2 X 10'^ 
(d)  An  = 6 X 10“^;(e)  pulse  incident  at  input  surface. 


index  variations.  Curve  e shows,  for  comparison,  the 
incident  pulse.  Curve  a shows  the  output  pulse  for 
An  = 0;  this  is  just  the  diffraction  spreading  that 
occurs  in  the  absence  of  a grating.  Curve  b,  c,  and  d 
show  the  output  for  increasing  values  of  An;  as  An 
increases,  the  pulse  becomes  increasingly  confined. 
However,  the  greatest  confinement  appears  to  be  of 
the  energy  that,  in  the  absence  of  the  grating  struc- 
ture, falls  in  the  middle  regions  of  the  spread-out 
pulse.  The  energy  in  the  outer  portions,  corresponding 
to  the  most  oblique  plane  wave  components,  shows  a 
lesser  tendency  to  confinement. 

The  plane  waves  satisfying  the  Bragg  condition  are 
represented  by  the  spatial  frequency  components 
±\jld,  which  is  only  4 of  the  maximum  spatial  fre- 
quency ltd.  Also,  we  note  that  the  value  of  An  re- 
ported by  Rosenberg  and  Chandross  in  their  experi- 
ments was  6.4  X 10'*;  this  is  close  to  the  value  of 
An  (6X10  ‘•f  for  which  our  curves  show  nearly 
complete  wave  confinement.  Finally,  we  observe 


that  thinner  structures,  such  as  dichrornated  gelatin, 
produce  refractive  index  changes  up  to  two  orders  of 
magnitude  greater  than  those  considered  here.  Thus, 
when  spatial  frequencies  on  the  order  of  a few  hun- 
dred cycles/mm  are  recorded,  the  guiding  conditions 
arc  generally  satisfied,  even  thouglr  the  guiding  effects 
arc  not  observable,  since  the  propagation  distances 
are  tor)  short  to  reveal  this  phenomenon. 

We  thank  the  Office  of  Naval  Research  for  their 
support  of  this  work  (contract  number  N 00014-67- 
A-0181-0058). 
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Pi  opagation  of  spatial  pulses  in  interferomeiricaliy 
produced  guiding  structures 

G.  C.  Kung  and  E.  N.  Leith 


Ir.iiKlcr  iir.'piTl  ics  iif  inlirlrnimct  rii  alls'  proilia  cil  diclii  l rir  ht;lil  uiii<li"s  arc  ms  csl  laal  ial  A L'l'iar 
a I I la  li  nil  pic  III  anals'sin  i^  ilrsi'lnpi'il  rmiipuli'r  laii'iilal  nm'i  liaaial  on  I Ins  a[i[irnat  li  an  nsial  In  si  ml'.  I Kc 
l||•llasnlr  III  (inlsis  irasi  rsini;  llic  pniilnii;  slrnilurcs  Inr  sariniis  salni's  nl  rcirail  is  c iinii'X  K.spirnni'iiial 
rosiili  s an-  nliiann  il  iisiiit:  plininsi-nsiii/cil  puls  mcllis  I iiii'l  Inn  rs  lain  as  i lie  rnrnriiiiii;  maicrial.  ( 'urn puli  r 
ri'sniis  arc  cinnparril  ssnii  l■xlstill(;  llinnry  anil  Innnil  to  lie  in  clnsn  at;rcfninnt. 


Introduction 

Optical  ss'iise  conliiicinpnl  in  dielectric  slalis  Imr- 
(lered  l)v  materials  of’  losser  rr-t'ract is'e  index  lias  lieen 
svidely  treated.  'I'he  general  treatment  ol  optical 
slat)  ss aseniiides  is  usually  liased  on  the  assumption 
that  tlie  cladding,  tlie  material  of  lower  relractise 
indt'x.  extends  to  infinity.  It  sescral  slali  svaveguitles 
are  placed  parallel  to  one  anothiT  in  space,  the  .1-- 
siimption  is  Usually  made  that  the  cladding  is  sulfi- 
eientiv  thick  that  the  presenc(‘  of  one  guide  does  not 
distort  the  field  of  otliers,  and  tliere  is  thus  110  inter 
action  lietsveen  the  ssas cgiiides.  'I'his  is  the  treat- 
ment int rodneed  liy  Tien.'  liossever.  in  mans'  prac- 
tital  lases,  tlic  assumption  tliat  parallel  slab  ssase- 
guides  ar(>  independent  of  one  anotlier  lannot  lie 
made;  such  svasegnides  must  lip  ( imsidered  to  lie  eon- 
pled.  'I’liis  is  the  case  svlien  the  slab  svas'cguides  are 
in  dost’  proximity  and  es))e(  i;iils'  svlien  lliere  exists  no 
disiri-te  lioundary  between  adjacent  svavegnides,  as 
oceiirs  svlien  they  are  formed  by  recording  ol  interfer- 
ence Iringes  in  a thick  dielectric  mediimi.  'Flic  con- 
[lied  svas'c  etpiations  lor  this  i ,-ise  het  onie  .so  iinman- 
agealily  s’oniplisatf'd  lliat  a ditferent  a)iproach  is  in 
order,  A rav  optics  approacli  was  used  liv  Hosen- 
herg  in  Ins  aii.ilvsis  of  llip  inlerleromel  riealls'  pro- 
diu'ed  Idler  opiics  slriu  tures  prodiii'cil  liv  Kosenlierg 
and  rii.'indross, 

l‘A  ideni  Iv.  an  arras  of  slali  svavegnides  formed  in 
tlii.-  Ill, inner  can  eipialis'  svell  lie  regard(-d  as  a solume 
pliase  grating,  and  we  sliould  exjieit  tlial  tlie  sarioiis 
anaivses  tliat  liave  lieen  given  lor  diltrac'tion  from 
volnnip  gratings  produced  by  recording  inlerterencc 

Till'  iiulliiir-  Iiri  ssilli  Ihr  I Iniscr-ils  ul  Mirliii;;iii,  Ann  \rlmr. 
-Mil  hiUMti  tsinT  • 

fi*.  I (I  M .Iiit.r  I !*7f» 


fringes  should  also  descriiie  image  transfer  effects 
manilesti'd  by  these  gratings.  We  aim  to  develop  a 
general  tiieory  tliat  s'ould.  liy  no  more  than  a change 
of  parameters,  desirilie  Bragg  diffraction,  tlie  image 
transfer  effects  that  occur  svilh  sufficiently  higli 
index  modnlatioii.  and  tlie  transition  ease  svliere. 
svilh  increasing  index  modulation,  the  Bragg  diffrac 
lion  elfeel  lieeomes  .ui  image  transl’er  or  guiding  el' 
lect. 

Method  ot  Analysis 

We  base  adapted  the  tliin-grating  decomposition 
method  of  .•Mieriiess  to  tlie  aiialssis  of  guiding  in  in 
lerleroniel  rieally  Inrnied  slali  waveguides  ''  Tliis 
nielliod  is  eonceptnally  ipiite  simple,  vet  it  gives  re- 
sults that  agree  exailly  svilh  tliose  o(  more  eomplex 
theories  suili  as  gis'eii  liy  Burckhardt,"  Kogeliiik,' 
atid  ixaspar.  ‘ I lie  melliod  consists  ol  dt*coni|iosing  a 
thick  grating  into  a -lacked  .irray  ol  I hi  11  gral  i ngs  and 
treating  eai  li  lliin  grating  somponent  m Ir  rnis  of  the 
svell  knosvn  tliin  pli,i-e  grating  tlKuirs'.  'I'lie  ditlrai  t 
ed  orders  ot  the  itli  thin  grating  tonus  tlie  incident 
svaves  ol  the  i + 1th  thin  grating.  Fhis  process  i- 
eonliniied  until  sve  n ,u  h tlie  linal  pi, me,  ss  here  tlie 
various  plane  ssas'p-  enu-rge  as  the  dilfraciid  order- 
ol  the  thii'k  grating  Tlic-e  svaves  are  just  tlie  otie- 
Ihat  salislv  tlie  Bragg  angle-,  altliougli  no  Bragg  dil- 
Iraclioii  con-ideral  ion-  sven  explicitly  in-erted  into 
tlie  llieorv.  I’lie  saiiinis  diffraeled  order-  that  do 
not  salisis  the  I'iragg  condition  simpis'  as'crage  to 
s .iliie-  near  /ero. 

\\  e liave  niodilied  the  .-Mli'niess  melliod  li\  using  a 
narriisv  di-l rilinl ion  of  light  (a  spatial  piilsel  as  the 
incident  Held  1 Fig.  11.  W e d(*eompo-e  I hi- pul-e  into 
.1  -pei'lrum  ol  plane  svases  and  propagate  e.uli  one 
indi'pendenlly  llirongh  tlie  -irueinre;  eai  li  svave  gen 
crate-  \ariou-  dlllr.u  led  orders.  All  svave-  are  then 
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I Pr''|M^:ntion  of  a pulse  of  li^hl  thmu^li  a Ihli  k ^'rafin^: 
».P  slunvs  the  rat  in;:  si  met  u re.  wit  h >.pa(  in^'  d:  i h > shows  the  input 
pulse,  anil  (r)  shi>ws  tlu*  pulse  emer^in*:  from  the  othi-r  surlare. 


summod  at  the  output  surface  to  yield  tlie  resultant 
intensity  distribution.  The  theory  as  used  here  in- 
volves the  same  assumptions  and  conditions  as  im- 
posed hy  .Mferness.  The  polari/ation  is  with  the  E 
vector  parallel  to  the  surface  of  the  recording:  materi- 
al. and  the  ditfracted  waves  are  summed  coherently 
at  each  successive  thin  crating:,  with  the  assumption 
of  no  polarization  changes,  followinw  the  earlier  treat- 
ments ol  Hurckhardt  and  Ko^felnik.  The  intensity  at 
the  exiiinc  plane  is  calcidated  without  regard  to 
propaiiat  ioti  direct  ion. 

■Mthouch  the  method  is  not  restricted  to  any  spe- 
cific prolile  of  input  pulse,  we  have  assumed  one  with 
amplitude  ol  the  form 

/Tt  ’•  ,n  ‘ c s2-  . ■■./)  rrci , V/,  (1 ) 

where  i!  is  the  cratinc  period  ilO  " mm  in  this 

р. iper),  and  is  ,i  const, int  that  adiusis  ihe  width  of 
the  input  pulse.  Not  only  is  this  pulse  spatially 
hounded,  its  spatial  frequeni  v spei  irum  is  .also  near- 
ly hounded.  For  compulation.il  puipo>es.  we  irun- 

с. ile  Ihe  spectrum  at  /,  = ±'2/'ui.  thus  ignorin';  the 
small  secondary  lobes  lyinc  outside  this  hound  I'his 
leads  to  only  very  small  errors  in  the  results  hec.uise 
the  secondarv  lobes  are  of  the  order  of  10  ’ that  of 
the  peak,  as  shown  in  Fic,  2.  In  our  an.ilvsis,  we  rep- 
resent /<v)  hy  a Fourier  series:  each  term  of  the  series 
represents  a (ilane  wave  that  we  can  pro[tan;ile 
through  the  f'ratim;-  I his  Fourier  series  represenia 
tion  results  in  afiditional  pulses  /'i.v  - nl‘),  where  /'  is 
the  period  of  the  lutidametftal  component  of  the  Fou- 
rier series.  We  must  choose  /’  siil  I i<  i<>nl  Iv  lar^e.  i.e., 
lake  m.'inv  samples  of  Ihe  spectrum  of /l  v ) so  that 
Ihe  various  input  pulsies  renuiin  suflicienlly  sepa- 
rated that  they  do  not  interact  as  thev  travel  lhrout;h 
I he  ural  Ilia. 

omputer  (iroeram  developed  hv  .Alt  •mess  for 
his  thin  (.'ratiiif;  decom[iosil ion  ;in.ilysis  was  modified 
for  tl>e  study  of  imat'e  transfer  efteds  The  comput- 
er profjram,  which  assumes  unsl.uitefl  friut;es  and  ;i 
sinusoidal  variation  of  refractive  index,  allows  for  ad- 
justment of  input  pulse  width  and  eraliiif;  thickness, 
later, d positioninu’  ol  the  initui  pulse,  selection  of  re- 
fractive index  modulation  n|.  anil  afijustmeni  ol  Ihe 
(,'ral infi  period.  'I'he  pulse  s[)ecl rum  is  compulerl  and  * 


sam[)led;  each  sample  represents  a [)lane  wave.  We 
can  select  the  numher  of  samples;  the  more  we 
choose,  Ihe  fjrealer  the  separation  of  .adiaceni  inci- 
dent pulses  iind  the  thickir  the  trratiuf;  that  can  he 
treated  vvithuvit  errors  Irom  interaction  helween  adja- 
cent pulses.  W e can  also  choose  th<*  numb,  r of  thin 
f;ratin(js  into  w Inch  the  thick  rat  in;'  is  .!••(  ■..•|,,isi  d 
We  follow  two  t'ener.il  f;uldehnes  in  i ho., >111;;  the 
numher  of  |ilane  wave  sam|)les:  the  thickness  of  the 

thin  t;ratin;;s  and  the  numher  of  diffracted  orders  to 
he  retained  from  e;n  h thin  gralin;;.  First,  the  [tres- 
ence  of  line  oscillation  patterns  at  the  ed;tes  of  Ihe 
s[)readoul  pulse  indicates  that  adjacent  [jiilses  are  in 
lerferini;  with  each  other,  and  theretore  more  plane 
wave  samples  are  needed.  The  heeinnine  of  such 
overlap  is  evident  in  Fi;;.  curve  a.  Second,  the 
total  enertty  at  the  input  and  out[)ul  planes  should  he 
conserved;  conservation  failure  indicates  that  thinner 
sections  are  needed  or  that  more  diffracted  orders 
should  he  kept.  'I'hroughout  our  analysis,  we  keep 
five  diffracted  orders  (0,  ±1.  i2  orders)  and  use  a 
thin  grating  thickness  of  tl.O.t  mm.  'I'hese  values  are 
generally  adetpiate  for  index  modulation  ft;  up  to 
0.001  ;■).  ’ 

Computed  Results 

Using  the  theory  thus  developed,  we  hiive  comput- 
ed the  behavior  of  pulses  traversing  the  guiding 
strui'lure.  Figure  0 shows,  for  a representative  con- 
liguration,  the  result  of  increasing  the  index  modula 


I' lU  < lnton-tl\  nl  output  surtaoc  lor  various \alurs 

ft  M r (a)  u 1 = ti.  I 111  M I = 10  *,  1 'j , = 2 X 10  ((i  1 •! , = (1  V 

10  tai  puUi'  mi  jdrnt  .11  input  surlme  (Iratuii:  thu  knrv*;  = 0 .M 
?nm 
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D I ST4\Cf  IN  MM 

Kit  I Samp  as  Kit.  't.  pxcp|iI  lor  a |)iiKp  wiiilh  of  0 0:!  rnni:  lal 

shows  Ihp  itu  ident  pulsp;  (h)  siiows  fhp  out|ni(  pulse  for  n \ = 0. 


PIS  T A N C f.  ! N M M 


Kit  a ( 'onl ioual ion  ot  Kit  !.  sliowint  oulput  pulsp  for  variims 
values  ol  ft  i‘  lat  ft  1 = 10  *;  |!>!  u l = :!  X 10  *;  (i  ) rj  i = l!  X 10 
Oil  0,  = 10 


lion.  As  the  index  increases,  tlie  pulse  width  at  the 
outffut  approaches  in  width  the  input  pulse,  and  tlie 
energy  in  adjacent  channels  diminishes.  In  P'ig.  li,  as 
well  as  in  succeeding  figures,  the  regions  of  maximum 
index  are  at  I),  -tO.OI  mm,  etc. 

Ke|)eating  the  same  experiment  with  incrr'ased 
pulse  width  produces  the  results  shown  in  Figs.  4 and 
The  pulse  width  is  now  O.O.T  mm,  so  that  a portion 
of  the  pulse  extends  over  the  houndaries  into  adja- 
cent channels.  As  n | increases,  the  output  pulse  is 
reshaped,  dividing  into  several  pulses  occupying  sev- 
er;d  adjacent  channels,  much  as  if  the  pulse,  like  an 
incom(>ressihle  fluid,  wen*  being  pressed  into  a mold, 
although  it  is  seen  that  the  reshaped  energy  dist rihu- 
lion  floes  not  exactly  conform  to  the  molfl.  We  note 
t hat  the  energy  dist  ri hut  ion  in  t hr-  various  channels  is 
not  alwavs  centered  at  the  ccnier  ot  the  channel, 
-^which  we  lake  to  he  the  points  of  maximum  index, 
occurring  at  0,  O.Pl,  (1.02,  . . . mm;  such  asvmmetrv  is 
related  to  tlie  asvnirnelry  of  the  impresserl  irradiance 
acrfiss  each  channel  input.  Finallv,  we  note  that,  de- 
spite the  over  all  confinement  of  eni'rgy  to  those 
channels  upon  which  (he  energy  was  incident,  there 
is.  with  increasing  fii,  an  increasing  amount  of  energy 
in  the  outer  channels,  whereas,  for  0|  = 0,  thi^Tc  is 
negligible  energy  in  this  region. 


Next  we  try  a pulsp  of  O.OO.'i-rnm  width,  which  is 
half  of  the  width  of  the  guiding  channels.  Figure  6 
shows  the  sequence  at  the  output  as  a function  of  rii. 
For  sufficiently  large  fi|,  the  |ndse  width  converges 
again,  as  in  Fig.  2,  to  that  of  the  input  pidse.  This 
result  is  representative  of  our  findings;  (or  a pulse 
narrower  than  a channel  width,  the  output  pulse 
width  converges,  with  increasing  fi|,  to  that  of  the 
input. 

If  the  pulse  is  introduced  at  a position  displaced 
from  the  center  of  a channel,  the  confinement  de- 
creases; the  losses  are  greatest  when  the  [ndse  is  cen- 
tered at  a |)oint  of  minimum  index  (Fig.  71.  The  dy- 
namics of  the  off-centered  pulse  are  further  shown  in 
Fig.  8,  which  traces  the  course  of  such  a pulse.  'I’he 
(lulse,  as  it  progresses  through  the  guiding  structure, 
oscillates  from  one  side  of  the  channel  to  th<“  other. 
As  the  pulse  swings  to  one  extreme,  some  energy  sur- 
ges across  the  refractive  index  harrier  into  the  adja- 
cent channel:  this  energy  produces  in  that  channel  a 
propagating,  laterally  oscillating  pulse,  which  in  turn 
impels  some  of  its  energy  into  the  next  adjacent 
channel. 


IN  lOM 


Kijj.  (!.  Samp  as  Ki|,'.  ■(,  pxd'pl  lor  a pulsp  width  of  0il0,'i  nm:  l.al 
ti  I = (I;  111)  n 1 = 10'  II  I It  I = g X 111  (d  l a , = ■(  x 10  (p|  a i = 

■H  X 10  (I)  a,  = g X 1(1-'. 
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C d 

l-'a;,  7 Klfpi  I ol  inIrodiH  iiii;  ihp  pulsp  al  an  oil  i-pnlprrd  posilion: 
piil-p  widlli  = 001,  ai  = 0 X 10  thiikiii'ss  - 0;',(  nim.  Dis 
plaiptiipiil  Irom  i cno  r i.U  0,  ihl  0 OO'J  mm,  li  ) o,oo:i  mm;  (dl 
0 00.'.  mm  \ prill  al  -i  alp  douhlvd  lor  0 1 and  id  I \'prliral  sralp  is 
inti'iisiix . hori/onlal  sc, dp  1 1 1 is  off  i rnli-r  dtspl.n  pim-nl  in  nun 
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Km  X Tracing:  nn  nl)  criilrr  injivH  piiUr  !hri>tu:h  llu* 
clunncl.  I h'^pl.m  nwtil  a!  im  idtMu  «•.  (i.lHl.i  mm.  pulse  viidlh  = 
0 nl  mm.  j = d X 10  ‘ /is  intensity.  » i'^  displ.n  ement  in  mm. 
Depth  at  (d)sprvnt  ion  plane  is  (a)  0.  (ht  tt  IH  mm.  let  0..'ir)  mm.  (cD 
O.al  mm. 


1' m 0.  Kner^jv  transfer  between  ( h.'uinels.  'I'he  lower  cur^e 
slu'ws  thr  pet<eni  decrease  ifi  ener^^y  in  \aricnis  ihannels  with 
propimalioti  distance-:  la)  is  lor  < linnnel  o itfu*  ch.innel  where  the 
pulse  was  introduced),  (b)  is  l«>r  the  + I channel:  and  It  ) is  tor  the 
-1  channel.  Shown  at  the  top.  for  reference,  is  the  deviatiiui  of 
the  pulse  from  the  channel  center 


lern  is  tlu*  same  as  lor  llio  (low  into  lire  ±1  cliamiels. 
and  the  previous  ol)ser\ations  apply.  For  compari- 
son, we  plot  (Fifj.  II)  the  enerfjy  <>ih  <>f  'lie  0th 
channel  for  the  case  where  the  pulse  is  introduced  at 
the  channel  center,  tlie  refractive  index  maximutn. 
As  expected,  the  eneri;v  loss  is  less  osrillalory,  al- 
thoufjh  a transient  occurs  at  llie  l)ej;intiint;.  Also,  the 
energy  leakage  is  significantly  lower. 

Experimental  Results 

Kxperimenlal  data  supporting  in  a (pialitative  way 
some  of  the  calctdated  results  were  obtained  hy  pre- 
paring polymethylmethacrylate  (I’MM)  casting  sen- 
sitized with  p-henzo  (piiiione:  at  present  we  do  not 
ha\e  sufficient  control  o\'cr  the  experimental  process 
to  prorl\ice  (junni it.il i\e  residts.  Using  the  ().-18,S-#jtn 
line  of  the  <»rgon  laser,  unslanted  fringes  of  lOO-l/mm 
sjiatial  frequency  were  reconled  in  O.o-min  thick 
specimens,  'f'o  demonstrate  image  transfer  charac- 
teristics. scratches  approximately  parallel  to  the 
fringes  were  made  on  one  face.  'I'he  PMM  was  then 
illuminated  with  a green  light  of  about  ,500  A in  spec- 
tral width,  with  the  light  iiniringing  on  the  scratched 


Km  10  Kn«T^;y  n»iw  int(»  nilu  r c}i.uincD.  slutwiim  oncr^y  iti  OiU'h 
( h.iniH-l  as  prrciuil  of  im  idriif  i ' (<d  is  ( banncl.  i b)  is  —2 

{■h.'innrl;  (cl  i"  +.1 1 luinm-l;  (d)  is  —'A  ch.tnncl. 


The  cur\es  of  Fig.  !)  show  that  the  greatest  rate  of 
energy  How  into  the  adjacent  channels  occurs  be- 
tween till'  positions  when  the  inilse  has  crossed  the 
ceiiler  of  the  channel  and  when  the  pulse  reaches  its 
maximum  olf-center  displacement.  'I'he  energy  flow 
into  the  two  channels  adjacent  to  the  0th  channel 
(the  one  where  the  pulse  was  introduced)  is  asym- 
nu  t ri^;al  Most  notalil>  , tlie  energy  liuildup  in  llie 
chaTmel  opposite  to  the  direction  o(  initial  pulse  dis- 
pl.'icement  (the  —1  channel)  ri'i'eives  a bigger  initial 
surge  of  energy  than  does  the  other  (the  1 channel), 
evidently  because  the  surge  of  energy  into  the  +1 
channel  is  aborted  hy  the  introduction  of  tlie  pulse  at 
a |)hase  (maximum  displacemen(  I jus(  pas(  (he  posi- 
(ion  of  maximum  energy  How.  Figure  10  shows  (he 
energy  How  into  other  channels.  .In  general,  the  |>at-  “ 


lOO 


0 d’  ' 40  60 


' -lie  “ \E  $S 

(■'n;  II  low'ruv  l,,v,i  whrn  ihi'  in|>ul  piiUi'  is  nil rmliu c<t  nl  tlie 
I biiniH-l  ccnicr  l.ii  vbnws.  flu  ju-ruf-nt  Itws  Imm  ch.mn<’l  (); 

lb)  «sb««Nvs  ilu'  in  tbr  + 1 or  - 1 r)i.inncl«i 
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Fi^.  12.  Oplicnl  for  oh.'ierving  pffecls.  0 is  the 

object,  a FM.M  ^Tatin^  with  scra(che.«  inscribed  on  the  left  surface. 
I.  is  a li^ht  beam,  ititroduced  oblit^uely  .so  that  no  (ilrect  liKht  prop- 
agates tiuough  the  system.  / is  the  image  of  the  PMM,  and  .^/  is  a 
microscope,  which  can  be  focused  onto  any  part  of  the  image. 


Fig.  n Kxperimental  results.  'Vhe  micri)scope  is  focused  on  a 
scratch  at  the  incident  surface.  The  upper  portion  of  the  photo- 
graph shows  the  scratch  in  sharp  focus;  in  this  regitm,  no  grating  is 
rec«*rded  on  the  FMM.  In  the  lower  portion,  the  scratch  is  out  of 
focus  due  tt»  the  presence  of  the  grating  waveguide.  In  the  center, 
the  miirosfope  was  focus(-d  on  the  exiting  surface;  here,  in  the  re- 
gion of  the  recorded  grating,  the  pulse  is  sharply  focused,  thus 
demonstrating  a strong  guiding  elfecl.  At  the  right  is  the  same 
image,  nvt*rexpo.''ed  so  as  to  emphasize  the  side-lolies.  i.e..  the  ener- 
gv  that  has  crept  into  other  channels. 


surl'ace.  This  illumination  was  infrodneed  ohlicnu'ly, 
so  that  only  scattered  li^jlit  entered  the  imaK'injj  sys- 
tem (Fifj.  12).  I'he  reeordiiiK  eamera  could  Ire  I'o- 
ensed  on  either  the  snuree  or  the  oliserver  side  of  the 
I’MM.  Results  are  shown  in  Fin.  11)  (which  is 
cro[)pcfi  to  show  only  (>ne  scratch).  At  the  left,  the 
.camera  was  focused  on  the  sonree  side.  The  record- 
ed fringes  are  confined  to  the  lower  portion;  thus,  at 
the  to)),  the  scratches  are  clearly  observed.  In  the 
renter  is  shown  the  same  imane.  I)ut  with  the  focus  on 
tlie  s\irfaee  toward  the  observer,  where  a stronn  fj'dd- 
ing  effect  is  evident.  'I’he  intensity  irregularities  re- 
sult from  nonuniformity  of  the  scratch.  At  the  right 
is  the  same  image,  using  greater  ex|)osure  on  the  re- 
cording film  so  as  to  eni|))iasizc  tlic  second.irv  or  sidc- 
lohc  structure  of  the  image.  I'he  lohing  region  is  of 
the  same  width  as  the  unfocused  beam,  with  the 
outer  lobes  Ireing  somewhat  more  intense  than  the 
inner  oni's. 

.‘since  the  scratch  is  not  aeenralely  ))arallel  to  the 
reeorrled  fringes,  there  is  a continuous  displacement 
as  a function  of  tlie  vertical  coordinate  he)^veen  the 
in|)ut  ))ulse  and  the  guiding  channels;  this  (auses  the 


|)ulse  ])eriudically  to  shift  channels.  'I'his  effect  is 
shown  in  more  detail  in  the  magnified  image  of  Fig. 
14.  Note  that  beginning  at  the  center  and  continu- 
ing downward  somewhat,  the  beam  evidently  falls  in 
a region  of  lowest  index  (i.e.,  straddles  the  channels), 
and  the  light  at  the  out|)ut  is  s|)lit  between  the  two 
adjacent  channels.  The  results  of  Fig.  1 t were  |)ro- 
duced  by  imaging,  using  a He  N'e  laser,  a narrow  slit 
onto  the  source-side  surface  of  the  RMM,  and  the  re- 
cording camera  focusing  on  the  surface  toward  the 
observer.  The  result  is  e.ssentially  no  different  from 
the  scratch  method  used  for  the  previous  results. 

Comparison  with  Previous  Results 

Rosenberg,  using  a ray  tracing  analysis,  has  com- 
|)uted  the  efficiency  of  image  transfer  in  interferome- 
trically  |)roduced  fiber  optics  structures.’’  Figure  4 
of  Ref.  1)  shows  guiding  efficiency  as  a function  of 
vertex  angles  of  incident  cones  of  light.  However, 
there  are  two  major  differences  between  his  case  and 
ours.  First,  we  treat  only  one-dimensional  guiding 
structures,  i.e.,  guiding  slabs,  whereas  Rosenberg  an- 
alyzes the  |)ro))erties  of  guiding  columns.  Second, 
our  incident  irradiance  is  a narrow  |)ulse,  having  most 
of  its  energy  concentrated  at  the  middle  of  the  guid- 


Kii:  1 I (bjitiinj;  nt  a slit  liuUsfd  onto  thr  FMM  svirlrtrr 

Iht*  sill  iv  sliijliily  iionp.ir.tlU'I  \Mtb  tho  jrrnlinc  frin^;<'s 

SonH'what  llu-  tnidtllr,  thr  puUo  i r»»ssrs  o\«‘r  thr  chnnnrl 

bouiul.irv  .uul  u split  inti'  tvvo  nr. iris  hju.i)  • Ttion« 
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Fij;.  IT)  Opt  ical  system  tor  ol>ser\ inn  (iil'fracliou  splitting,  / 1 ami 
I t are  the  liroadened  first  order  diffrai  led  inianes.  I‘,  is  the  soorie 
plane.  Source  is  imaned  onttJ  anti  finally  onto  I* t. 


iiig  channel,  where  (he  fritidinu  elliriencv  i.s;  hijjhest, 
whereas  Rosenberg  considers  rays  trf  etjual  intensity 
in  various  directittns.  We  define  guiding  efficiency 
as  the  percent  ttf  incident  light  that  remains  in  the 
0th  channel  after  propagating  through  the  guiding 
structure,  f'or  comparison  of  our  results  with  liosen- 
herg's,  we  choose  parameters  idetitical  to  those  of 
Rosenberg;  we  choose  Sn  (=  2rii)  = O.OOOti,  thickness 
of  the  guiding  structure  = 1.1  tnm.  and  fritige  spatial 
frequency  = 100  I/mm. 

From  our  computer  result,  we  find  th.it  the  guiding 
efficiency  of  a 0.01-mm  wide  pulse  in  the  form  of  Fa). 
(1)  is  42%.  Since  the  spectral  distribution  of  our 
pulse  is  not  uniform,  we  cannot  directly  compare  our 
result  to  that  of  Rosenberg.  From  the  s|)ectruin  of 
our  input  pulse,  we  determine  that  12‘Vi  of  the  total 
energy  lies  between  a wedge  of  ()..'1°,  I 1. more  lies 
within  0.fi°,  11.2'H>  more  within  0.0°,  lO.fi'y,  more 
within  1.2°,  and  so  on.  From  the  guiding  efficiency 
curve  of  Rosenberg,  we  can  estimate  the  guiding  effi- 
ciency for  the  in[)ut  irradiance  lying  within  the  vari- 
ous wedges.  Multiplying  these  efficiencies  by  the  per- 
cent of  total  input  energy,  we  estimate  the  over  all 
guiding  efficiency  to  he  l.V.’t.,  which  is  in  close  agree- 
ment with  the  result  predicted  by  our  analysis. 

Diffracted  Order  Splitting 

Another  verification  of  the  validity  of  the  thin 
grating  decomposition  analysis  is  its  success  in  pre- 
dicting a phenomenon  that  is  the  optic.d  countertiarl 
of  a well-known  x-ray  diffraction  effect  predicted  by 
Kato'”;  namely,  a plane  w.ive  incident  at  or  near  the 
Bragg  angle  of  a sinusoidally  stratified  structure 
splits  into  four  plane  waves.  We  report  some  addi- 
tional observations.  Instead  of  using  a narrow  beam 
incident  at  the  Bragg  angle,  as  Forshaw"  did.  we  illu- 
minated the  RMM  sani|)le  by  a s|),itiall\  incoherent 
slit  source  with  a beam  divergence  such  as  to  include 
rays  satisfying  the  Bragg  condition  (or  both  first  or- 
ders. The  slit  image  was  formed  about  I mm  in  front 
ot  the  I’MM  surfaie.  U'e  then  made  observations  on 
the  sLTflnage 

The  experimental  setup  is  shown  in  Fig.  I.'i.  The 
slit  source  was  prorluced  by  a Me  Ne  laser  beam 
passing  through  a rotating  grouiui  glass  placed  in 
front  of  a slit  aligned  parallel  to  the  grating  fringes. 
A unity  magnification  imaging  systi-m  consisting  of  a 
lens  pair  7^,  1-2  was  used  to  image  the  slit  into  the  vi- 
cinity of  the  I’MM.  A second  lens  pair  L\,l,.x  reim- 


aged the  P.MM  and  slit  image,  and  this  final  image 
was  observed  with  a mii  roscope,  which  could  be  fo- 
cused on  any  [ilane  in  the  vicinity  of  the  PMM  and 
slit  images. 

We  illuminated  first  an  area  of  the  PMM  sample 
on  w'hich  interferem c fringes  of  low  index  modula- 
tion were  recorded  .ind  uliich  produced  we.ik  dilTr.ic 
tion.  We  observed  th.it  the  first  order  dillracted 
beams,  as  observed  in  the  plane  of  the  slit  image,  ap- 
(icared  as  t wo  uniform  b;inds  of  light . As  we  changed 
the  distance  between  the  PMM  plate  and  the  slit 
image,  while  keeping  the  microsco()e  focused  on  the 
slit  image,  tlie  form  ol  the  light  distribution  remained 
the  same.  exce|)t  for  the  increase  in  distance  between 
the  diffracted  orders.  For  small  source-plate  dis- 
tances, th<‘  orders  overlapped;  our  observations  here 
apply  to  distances  sulficiently  large  that  the  orders 
are  distinctly  separated. 

A sini|)le  model  that  explains  the  observed  elfects 
is  shown  in  Fig.  111.  The  rays  that  satisfy  the  Bragg 
angle  produce  diffraction  as  they  propagate  through 
the  plate.  Rays  diffracted  at  different  planes  in 
depth  prodvice  diffraction  with  difterent  displace- 
ments from  the  zero  order.  For  the  weak  dilTraction 
exhibited  by  this  exposure,  the  zero  order  is  negligi- 
bly attenuated  as  it  propagates  through  the  plate, 
and  subsequent  diffraction  produced  by  the  diffract- 
ed waves  is  also  negligible.  These  assumptions  imply 
a uniform  intensity  distribution  across  the  diffracted 
order,  just  as  is  observed.  The  width  and  position  of 
the  diffracted  orders  were  measured  and  found  to  he 
just  as  predicted  by  this  simple  model. 


Kiu  Hi  Mntlrl  f<»r  tiImmi:  hm.iilt  nin>:  i»l  tiillr.it  li-ri  nnirr's 
U\  .Ttitl  A',  .irr  r.i\*»  at  llu  Hr.iui;  aiu(  /.•  :iri'  Uu  iir>st 

riillratti'd  ttrdcr'i  > is  lln‘  |>I.hh'  of  lh<‘  snircr  is  ihe 
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Ki^  IT.  I distrj!)uti(»n  Hi  iho  plaiu' ctl  Frirnnholer  ditfrMc  tuin 
(/'.  of  Hu'.  l'»).  showing  zerc»  orfl^'r  (li-fl)  ntui  om*  first  <>r(l»*r. 
('urvt's  a to  ft  reprfsont  incroasinfi  in<l*-x  tiMKltilation;  c is  a roin 
put er  lalfulnt ion  for  n i = 0 ot)()7. 


Illuminnt ing  rt'fjion.';  of  the  PMM  havint,'  larger 
index  modulations,  we  observed  that  the  images  of 
the  dilfraeted  orders  developed  |>eaks  at  the  edges 
(Fig.  17);  also,  a hand  of  light  appeared  in  the  vicinity 
of  the  zero  orders.  'I'hese  are  the  effects  predicted  hy 
Kato.  'I'he  latter  can  he  explained  as  being  produced 
hy  a sei'ondary  diffraction,  whereby  the  light  in  the 
first  order  is  rediffracted  hack  into  the  zero  order. 
The  migration  of  first  order  light  to  the  edges  of  the 
pattern  is  less  easily  explained  in  terms  of  an  intui- 
tive approach. 

'I’he  computer  analysis  |)redicts  exactly  these  re- 
sults. The  com|)uter  results  also  predict  the  width 
and  position  of  the  difiracted  orders.  However,  the 
coTjiputer  predicts  symmetry  about  the  center  ol  the 
grating,  whereas  exjterimental  results  show  ati  asym- 
metry that  we  presume  is  due  to  the  absorption  in 
the  I’MM;  the  absorption  loss  is  about  .'!()"<>.  We  find 
it  (pnte  interesting  that  an  anaivsis  based  on  nothing 
more  than  a continuing  itertition  of  the  (ptite  simple 
thin  grating  diffraction  theory  should  so  readily  pre- 
dict these  rather  subtle  effects. 

It  we  observe  planes  other  th.in  the  source  plane, 
we  find  that  the  spreadout  zero  order  j-ontains  two 
dark  bands;  these  corres[)ond  to  propagation  direc- 
tions that  satisfv  the  Bragg  condilion.  and  the  miss- 
ing light  is  thus  that  which  forms  the  first  orders. 
.Also,  :is  we  move  from  the  source  plane,  the  lirst 
order  diffracted  waves  become  bro.ader  and  the  edges 
less  shar|>;  this  is  expected,  since  it  is  an  angular 
banfl  ol  light  that  satislies  the  Bragg  condition,  the 
angle  being  determined  by  the  Bragg  angular  selec- 
I ivity,  w hich  depends  on  t lie  gral  ing  t hickness 

Simplified  Approach 

We  suggest  that  approximate  and  uselnj,  results 
can  be  obtained  for  the  sinusoidally  profiled  array  ol 


guifles  by  comparing  it,  on  a tpiite  simplified  basis, 
with  the  volume  grating  of  sinusoidal  [irofile.  De- 
spite the  sinusoidal  structure,  we  assume  that  in  the 
Bragg  diffraction  process,  all  rays  arc-  reflected  from 
discrete  locations  in  the  fringe  contours;  we  take 
thc'se  locaticais  to  be  at  or  near  the  n gions  of  mini- 
mum refractive  index;  we  assume  one  rellection  per 
fringe  period.  A ray  entering  the  material  is  thus 
partially  retlected  from  each  recorded  fringe  (Bragg 
surface)  it  encounters.  .Such  a simplification  has 
precedent,  being  the  usual  way  of  describing  the 
I.ippman  method  of  color  iihotography  and  a com- 
mon way  of  describing  holograms  made  by  Den- 
isyuk's  method. 

Confinement  of  rays  between  adjacent  fringes  oc- 
curs when  an  incident  ray  is  entirely  reflected  from 
the  first  Bragg  surface  encountered.  The  depth  of 
index  modulation  n\  recptired  to  produce  this  condi- 
tion can  be  obtained  from  the  relation 

/ C]  ciisib  2 t2) 

as  given  hy  Kogelnik,’'  where  is  the  grating  thick- 
ness, \|)  is  the  free  spaie  wavelength  of  the  light,  0|  is 
the  magnitude  of  index  variation  ttlui  = n,„;„  — 
n,,,,„).  and  ll,i  is  the  angle  of  the  incident  beam,  mea- 
sured within  the  material,  assuming  the  Bragg  condi- 
tion is  satisfied.  When  Kq.  (2)  is  satisfied,  the  dif- 
fraction efficiency  is  lOO",,,  that  is,  all  the  incident 
light  appears  in  the  diffracted  wave.  We  consider  a 
grating  of  thickness  7],  such  that  eat  h ray  encounters 
only  one  Bragg  surface  in  traversing  the  material; 
from  simple  geometry  we  find  that 

(ri;ui'>„  % t’l,.  (3) 

Combining  this  relation  with  the  Bragg  diffraction 
law  2d  sinfle  = Xd/uu  yields 

/'„  2l/f,r/'  A...  (4' 

where  n,t  is  the  bulk  index  of  the  material.  Fur  guid- 
ing, we  require  lOtr'.i  diffraction  efficiency  for  the 
thickness  7',,.  (hat  is,  a single  Bragg  reflection  jilacing 
all  the  incident  light  into  the  diffracted  order.  Thus, 
we  Combine  F,qs.  (2)  and  (4)  and  taking  costt,,  = 1,  we 
obtain 

II]  ft  4l/,)t.X(i  tl)~  n e tfi) 

as  the  minimum  refractive  index  modulation  re- 

(piired  for  guiding,  where  the  third  member  of  the 
eipiation  is  readily  obtained  from  thf>  s(“cond  by  ap- 
plication of  the  Bragg  ditfraction  law  The  reason- 
ableness of  t his  result  is  di'inonstrated  b\-  com|)arison 
with  the  comparable  condition  for  guiding  in  slabs 
with  sharp  boundaries.  ,At  the  boundary  between  an 
index  Ol,  — and  a higher  index  o,^  + n^.  total  re- 
flei  lion  occurs  lor  the  condition 

i‘i  ' i"  I'b’.  (6t 

a relation  readily  derived  from  Snell’s  law  and  a 
small  angle  a|)proximat ion.  I’he  two  relations  are 
nearlv  identic. d;  the  (actor  4 between  them  suggests 
that  the  reipiired  index  change  is  greater  for  sinusoi- 
dal va rial  ions  t h.in  for  discrete  boundaries, 

.Substituting  into  I'lq  t.'i)  the  values  assumed  in  the 
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computed  curves  (no  = 1.49,  Xo  = 0.6.128  ^m,  and  d = 
0.01  mm),  we  obtain  n\  = 6.7  X 10“'*,  which  is  indeed 
about  tl  e value  for  which  substantial  wave  confine- 
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Nonlinear  holographic  waveguide  coupler 


B.  J.  Chang  and  S.  K.  Case 


We  have  shown  ex[M*rimentally  that  tor  high  exp<isure,  higher  order  Hragg  diffraction  by  volume  holograms 
Is  dominated  by  the  nonlinear  refractive  index  modulation  and  proposed  the  a[)plication  of  nonlinear  modu- 
lation to  the  construction  of  a high  frequencv  volume-grating  waveguide  cou|)ler. 


Introduction 

[t  is  well  known  that  first -order  HraKg  reconstruction 
of  volume  phase  holoj(ram  gratinfis  (in  which  the  di- 
electric constant  is  spatially  modtilated  in  accordance 
with  the  exposure  intensity  distrihittion)  can  (jive  dif- 
fraction efficiencies  approaching  UKHi.'  Higher  order 
Bragg  reconstruction  of  volume  holograms  is  a lesser- 
ktiown  phenomenon,  however.  We  have  obtained 
strong  second  order  diffraction,  comparable  in  efficiency 
to  the  first  order  diffraction,  from  dichromated  gelatin 
holograms.  In  this  paper,  we  briefly  discuss  tbe 
mechanism  of  higher  order  Bragg  reconstruction  and 
propose  a unique  volume  phase  grating  coupler  whose 
construction  is  simplified  over  that  of  other  grating 
couplers. 

Higher  Order  Bragg  Diffraction 

Rose  and  Friesem  observed  strong  higher  order  Bragg 
diffraction  from  dichromated  gelatin  holograms.-  We 
found  experimentally  that  for  high  exposure,  higher 
order  Bragg  reconstruction  is  dominated  by  nonlinear 
modulation  of  the  refractive  index,  an  observation  that 
is  in  itgreemeni  with  the  [)redictions  made  by  Su  and 
(lavlord  ’ in  their  extension  of  coupled  wave  theory  to 
nonsinusoidal  gratings  and  to  higher  order  Bragg  angles. 

Higher  order  Bragg  diffraction  c:tn  arise  from  two 
sources:  first  order  diffraction  from  harmonic  gratings 
created  by  the  nonlinear  response  of  the  recording 
niitterial  to  the  exposure  interferenct  pattern  and  higher 
order  diffraction  from  the  ftmdamentjil  tind  other  lower 
orfler  gratings. 
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Although  higher  order  Bragg  diffraction  is  therefore 
the  summation  of  all  po.ssible  diffracted  orders  propa- 
gating at  the  appropriate  angle,  in  the  high  exposure 
range  the  contribution  of  lower  order  gratings  is  small; 
thus  the  Bragg  diffraction  efficiency  of  the  hth  order 
can  be  approximated  by-* 

’ll,  ~ sin'^lirnijd/tX,  cosW(,)|.  (1) 

where  X,  is  the  reconstruction  wavelength,  d is  the 
thickness  of  the  hologram,  m,  is  the  refractive  index 
modulation  of  the  /?th  harmonic  grating  when  the  re- 
fractive index  profile  ol  the  hologram  can  be  expressed 
as  either  a Vourier  sine  series  or  a Fourier  cosine  series, 
and  (I/,  is  the  Bragg  angle  of  the  hth  order,  determined 
by 

///.  - sin  , l'),  (21 

where  /i  is  the  fundamental  spatial  frequency  of  the 
hologram  grating.  F.quation  ( 1 1 is  an  accurate  ap- 
proximation under  certain  conditions';  otherwise,  this 
is  only  a rough  a()proximation.  .A  detailed  discussion 
of  the  mechanism  responsible  for  higher  order  Bragg 
diffraction  is  beyond  the  scojie  of  this  paper. 

'I'o  verify  experimentally  that  higher  order  Bragg 
diffraction  is  dominated  by  harmonic  grating  diffraction 
we  recorded  plane  wave  holograms  in  dichromated 
gelatin  derived  from  b49F  [ilates.  By  slightly  modifying 
the  processing  procedures.'  we  were  able  to  vary  the 
amount  of  nonlinearity  in  the  hologram.  If  higher  order 
diffraction  by  the  fundamental  grating  is  dominant,  we 
Would  expect  an  approximately  constant  second  order 
diffraction  efficiency  from  holograms  with  different 
amounts  of  nonlinearitv  and  the  same  first-order  dif- 
fraction efficiency.  We  observed,  however,  that  under 
these  conditions,  second-order  diffraction  efficiency 
varies  over  a wide  range.  In  fact,  with  the  appropriate 
nonlinear  modulation,  volume  phase  holograms  can 
diffract  nearly  all  the  incident  light  into  any  desired 
order. 

The  angular  bandwidth  of  the  hth  diffracted  order 
of  a volume  grating  having  a sinusoidal  index  modula- 
tion is  approximately  proixirtional  to  1/h  times  the  ratio 
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of  the  fringe  spacing  and  the  hologram  thickness;  this 
is  evident  from  an  examination  of  the  coupled-wave 
theory  and  /; -vector  diagram.'  Therefore  the  angular 
bandwidth  of  the  htb  (Mffracted  order  is  approximately 
the  same  whether  it  is  the  first-order  diffraction  from 
the  hth  harmonic  grating  or  the  corresponding  higher 
order  diffraction  from  the  fundamental  and  (»ther  lower 
harmonic  gratings,  i.e.,  the  angular  bandwidth  of  hth 
diffracted  order  is  apprcvximately  1/h  times  that  of  the 
first  order  of  the  fundamental  grating. 

In  Fig.  1,  we  show  diffraction  efficiency  curves  of  the 
first  and  second  orders  reconstructed  at  their  appro- 
priate Bragg  angles;  the  holograms  were  recorded  at  a 
wavelength  of  0.514  am  and  monitored  at  two  wave- 
lengths. The  diffraction  efficiency  of  the  second  order 
resulting  from  a large  refractive  index  modulation  of  the 
harmonic  grating  can  approach  that  of  the  first  order; 
indeed,  we  have  measured  efficiencies  that,  corrected 
for  reflection  loss,  are  in  excess  of  9,5"«i. 

Figure  2 shows  experimental  curves  indicating  the 
angular  bandwidth  of  the  first  diffracted  order  at  a re- 
construction wavelength  of  1 .06  ^m  and  of  the  second 
diffracted  order  at  a wavelength  of  0.514  am.  I'hese 
diffracted  orders  have  approximately  the  same  Bragg 
angle.  The  angular  bandwidth  of  the  first  diffracted 
order  at  a wavelength  of  0.514  am  is  of  course  similar  to 
that  ofthe  first  diffracted  order  at  a wavelength  of  1.06 
am,  but  the  Bragg  angles  are  different  (assuming  a 
constant  ratio  of  refraction  index  modulation  /i|  and 
readout  wavelength  X,  ).  In  Fig.  2,  we  see  that  the  sec- 
ond-order curve  is  approximately  half  of  the  width  of 
the  first-order  curve. 

Nonlinear  Grating  Coupler 

Kogelnik  and  Sosnowski  describe  a volume  holo- 
graphic coupler  as  a convenient  means  to  feed  light  into 
a thin  optical  film.'’  For  the  construction  of  this  volume 
holographic  coupler,  they  use  a prism  couitler  to  intro- 
duce the  construction  beams  at  an  angle  greater  than 
the  critical  angle  for  the  waveguide.  For  .some  materials 
and  devices,  it  may  he  undesirable  to  place  a prism  and 
index  matching  fluid  in  contact  with  either  the  photo- 
sensitive layer  or  the  substrate.  We  alleviate  this 
problem  by  exposing  the  film  with  one  pair  of  external 
plane  waves  and  using  a nonlinear  film  response  to 
generate  a harmonic  grat  ing  for  which  one  of  t he  read- 
out waves  will  he  a guided  wave. 

;\s  shown  in  Fig.  2.  ex|)osure  with  plane  waves  .4  and 
fi  of  wavelengt  h X forms  a grat  ing  characterized  by  the 
grating  vector  K|.  which  is  perpendicular  to  the  fritige 
planes  and  has  a magnitude  K[  = 2;r/  A.  where  A is  the 
f ringe  spacing  of  the  fundamental  grat  ing.  The  Bragg 
e(|uat.ion  of  the  /ith  ditfracted  order  is  then  given  hv 

sinB,,  = K'„  = hK  C(i 

where  ,i  = 2n-n,i/X  is  the  |)ropagation  <-onstant  in  the 
recording  medium,  (h,  given  hv  F,q.  (2)  is  the  angle  in  the 
recording  medium,  and  n,,  is  the  hulk  refractive  index 
of  t he  recording  medium. 

If  the  filn)  has  a nonlitiear  response  to  the  incident 
illumination  at  high  exposure  levels,  a second  order 
grat  ing  K . = 2K  i of  t w ice  t he  f undamental  sp.it  iai  f re- 
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qiiency  will  be  formed  within  the  emulsion.  I he  a- 
veetor  diagram  in  Fig.  4,  for  correct  construction  ge- 
ometry and  readout  with  wave  A',  indicates  how  we 
produce  a Bragg  diffracted  wave  H'  at  an  angle  greater 
than  the  critical  angle  for  the  guide  (i.e.,  a trapped 
wave).  Since  the  gelatin  and  the  glass  substrate  have 
approximately  the  same  refractive  index,  the  trapped 
wave  is  guided  in  both  the  film  and  the  substrate.  Since 
the  substrate  thickness  t » X,  our  guide  is  multimode. 

In  order  to  verify  the  above  theory,  we  constructed 
a grating  coupler  as  shown  in  Fig.  .I  with  plane  waves 
incident  at  and  70°  external  angles  with  respect  to 
the  film  normal.  The  polarization  of  light  for  both 
construction  and  readout  is  perpendicular  to  the  plane 
of  incidence.  The  readout  geometry  of  this  grating  is 
shown  in  Fig.  5(a).  A collimated  beam  of  coherent  light 
incident  from  the  upper  right  passes  through  a O.ITcm 
X 2-cm  aperture  and  illuminates  the  1 -cm  X I -cm  cou- 
pler grating.  I.ight  diffracted  by  the  grating  is  removed 
from  the  incident  beam  while  the  undiffracted  light  falls 
on  a piece  of  ground  glass.  Figure  5(b)  shows  the  sys- 
tem with  the  coupler  rotated  away  from  the  Bragg  angle, 
so  that  there  is  no  diffraction.  Figure  5(c)  shows  the 
grating  at  its  .second -order  Bragg  angle  and  nearly  100% 
coupling  efficiency,  as  is  evident  by  the  depletion  of  the 
incident  wave  in  the  coupler  region.  In  this  figure, 
however,  the  guided  wave  again  strikes  the  coupler 
grating  after  its  first  guided  bounce  and.  since  it  is  at  the 
complementary  Bragg  angle  for  the  coupler  grating,  is 
uncoupled  with  nearly  lOO'V  efficiency.  If  the  grating 
is  translated  so  that  the  guided  wave  will  miss  the  cou- 
I)ler  grating  after  its  first  bounce,  the  light  can  be  seen 
in  this  time  exposure  because  the  gelatin  of  the  film 
■ catters  a small  amount  of  light.  Removal  of  the  gelatin 
from  areas  of  the  guide  where  there  is  no  coupler  elim- 
inates this  scatter  loss. 

This  coiiftler,  of  course,  could  also  be  used  in  con- 
junction with  a thin  film  guide  for  single  mode  opera- 
tion, but  by  having  the  thick  guide,  we  have  gained 
coupling  efficiency*’  and  have  greatly  relaxed  the  re- 
quirement on  the  sharpness  of  the  grating  coupler  edge. 

The  authors  th.ink  K.  \.  Keith  and  VV.  S.  ('olburn  for 
hel()ful  discussions  during  the  course  of  this  work.  We 
thank  the  Air  Force  Avionics  Laboratory,  the  Office  of 
Naval  Research,  and  the  National  Science  Foundation 
(grant  OK  151  IH)  for  their  generous  support. 
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(*rtur(*  A and  illuminates  the  nmpler  ^ratln^  which  is  supported 
hy  a 1 - mm  thick  j'lass  substrate  l\  'The  transmitted  li^ht  illuminates 
^iround  ijlass  t h)  Illumination  tar  from  the  Hra^^  an^rle  ol  <-oupler 
^rat int:  let  Brazil  an^le  illuminat ii*n  near  the  center  tit  the  coupler 
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1 . Introduction 

Present  couplers  for  integrated  optical  circuits  ( 1 ] 
have  been  restricted  to  coupling  one  external  wave 
into  one  guided  wave.  We  concern  ourselves  in  this  re- 
port with  a thick  holographic  grating  coupler  designed 
to  provide  simultaneous  excitation  of  more  than  one 
guided  mode  when  illuminated  by  a single  external 
beam.  Once  the  guided  waves  have  been  processed 
within  the  guide,  we  will  want  to  selectively  decouple 
one  of  the  guided  modes  without  perturbing  the  other. 
A single  grating  decoupler  is  designed  for  this  purpose. 


2.  Double-grating  coupler 

We  have  shown  earlier  [2]  that  it  is  possible  to 
double  expose  a thick  grating  such  that  it  acts  like  a 
double  beam  splitter  to  diffract  50'/t,  of  the  incident 
energy  into  each  of  two  diffracted  orders.  We  use  this 


t le.  t . Consttuclu'n  of  a double  (iraime  coupler  A.  B.  ( ate 
mcidenl  ptane  waves;  1-  is  ihe  dichroinaled  gclalm  emulsion 
on  substrate  S,  P is  an  absorptive  backinn  plate.  Xylene  is 
used  as  an  index  matehint-  lluid  between  Ihe  componenis. 


principle  to  construct  a two-mode  waveguide  coupler 
as  in  fig.  1 . A dichromated  gelatin  film  is  exposed 
through  a prism  first  with  plane  waves  A and  B,  and 
then  with  plane  waves  A and  C.  After  development, 
readout  with  common  wave  A will  produce  both 
guided  waves  B’  and  C (fig.  2)  so  that  we  can  simul- 
taneously excite  two  modes  in  the  guide  with  a single 
incident  beam. 

We  have  also  constructed  two  single-grating  de- 
couplers by  exposing  a piece  of  film  through  the 
prism  with  beams  A and  B,  and  by  exposing  another 
piece  of  film  with  beams  A and  C.  Each  single-grating 
decoupler  is  exposed  and  developed  on  its  own  1 mm 
thick  ^ass  substrate  and  can  now  be  mounted  on  a 
waveguide  to  couple  or  decouple  beams  of  light.  By 
choosing  the  construction  waves  B and  C sufficiently 
far  apart  (in  angle),  we  can  construct  the  single  grating 
decouplers  such  that  they  act  independently  on  either 
of  the  two  guided  waves  13). 

In  fig.  3a,  we  show  a double-grating  coupler  mount- 
ed on  top  of  a i"  Kodak  Microflat  glass  guide  (emul- 
sion removed)  using  xylene  as  an  index  matching 


l it-.  2.  Rcadoul  of  llu-  doublc-)itatint!  coupler  with  wave  A' 

produces  guided  wavc.s  B and  C in  guide  G. 
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I 1^'  3.  I><»uble*^ra(ini:  coupler:  a)  Two  guided  modes  are 
simultaneously  excited  in  the  wave>!uide.  b)  One  ot  tlie 
euided  waves  is  decoupled  by  a sin^le-eratinp  decoupler.  i ) 
Moth  ot  the  I'uuled  waves  are  decjuipled  by  separate  sinele- 
>:ratint!  decouplers. 

(liiiil.  The  col!inut(t\l  beam  of  ligitt  (boaiti  A' ) incident 
trnm  Ilie  upper  right  is  coupled  with  iicrirly  100';  ef- 
I'icieticy  into  the  guided  beams  B'  and  C'  as  observed 
by  the  lack  of  transmitted  liglit  on  tlie  piece  ofground 
glass  below  the  guide,  Ilte  two  guided  modes  can  be 
seen  in  this  view  after  they  have  been  uncoupled  by  a 
prism  and  illuminate  a piece  ofground  glass  held  at 
grazing  incidence.  In  fig.  ,fb,  we  ailtich  a single-grating 


decoupler  to  the  bottom  of  the  guide.  This  grating 
was  constructed  with  waves  A and  C so  that  on  read- 
out. it  will  decouple  guided  wave  C'  and  difTact  it 
into  the  direction  of  wave  A'.  Guided  beam  B’  is  not 
at  the  Bragg  angle  for  this  decoupler  and  remains 
trapped.  In  fig.  .^c,  we  place  the  single  grating  con- 
structed with  waves  A and  B on  top  of  the  guide. 

Tliis  now  decouples  the  energy  from  the  remaining 
guided  wave  B'  into  the  unguided  wave  A’.  The  out- 
put beams  can  be  diffracted  in  either  an  upward  or 
downward  direction  by  a suitable  rotation  of  the  de- 
coupler gratings. 

3.  E.xtensions  of  the  principle 

Tliis  method  of  double-grating  coupler  construc- 
tion (with  common  reference  wave  A)  requires  that 
the  construction  and  readout  wavelengths  be  the  same 
so  that  the  two  gratings  will  have  a common  Bragg 
angle  on  readout.  If  the  coupler  is  to  be  used  at  a 
wavelength  different  than  that  of  construction,  then 
the  geometry  of  the  construction  beams  must  be  ad- 
justed so  that  the  two  gratings  will  have  a common 
Bragg  angle  at  the  readout  wavelength. 

By  making  the  mode-selective  decoupler  gratings 
sufficiently  thick,  their  angular  selectivity  can  be  re- 
luced  so  that  they  extract  light  over  a very  narrow 
angular  range.  In  this  way.  the  multi-mode  coupler 
and  decoupler  system  can  be  extended  to  include 
many  more  guided  waves. 
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We  advance  here  the  concept  of  modular  optics  and  its 
application  to  waveguide  optical  systems.  We  describe  a 
modular  optical  element  as  an  element  constructed  on  ins  own 
small,  thin  substrate'  and  capable  of  being  selectively  attached 
to  the  surface  of  an  optical  waveguide  and  later  removed. 
These  optical  elements  can  be  mass  produced  and  individually 
tested  before  use  in  an  optical  circuit.  In  addition,  if  the  el- 
ements are  holographically  formed,  they  do  not  have  to  be 
developed  while  attached  to  the  waveguide,  thus  saving  the 
guide  from  possible  chemical  attack. 

We  previously  used  thick  hologram  optical  modules  to 
couple  and  decouple  waves  selectively  from  a thick  wave- 
guide.'^ After  a coupler  module  was  attached  to  the  guide  to 
excite  two  guided  waves  within  the  waveguide,  thick  emulsion 
modular  decoupler  gratings  were  selectively  attached  to  the 
guide  at  positions  where  we  wished  to  sample  one  of  the  guided 
waves.  In  a similar  manner,  Wei  and  Tan  * have  developed 
a thin  plastic  grating  coupler  that  can  be  fixed  to  a guide  and 
later  removed. 

In  Fig.  1,  we  show  a construction  geometrv  for  a holoeraphic 
optical  module.  Wave  A in  Fig.  1 is  nearly  normal  to  the 
emulsion  E of  the  module.  Waves  B and  C are  at  angles  such 
that  they  will  produce  guided  waves.  We  first  construct  a 
coupler  grating  by  exposing  a film  module  with  waves  A and 
B.  Readout  of  this  mrxlule  with  wave  A will  produce  a guided 
wave  B.  We  also  construct  a mode  converter  by  exposing 
another  piece  of  film  with  waves  B and  C.  When  placed  on 
a guide,  this  module  will  allow  diffractive  coupling  between 
waves  B and  ('. 

After  development,  the  coupler  module  is  mounted  on  top 
of  a planar  gla.ss  guide  using  xylene  as  an  index  matching  fluid. 
In  Fig.  iha)  the  beam  incident  from  the  upper  left  is  diffracted 
into  guided  wave  B by  our  coupler  module.  Wave  B is  later 
decoupled  by  a prism  and  strikes  a card  held  at  grazing  inci- 
dence. 

In  Fig.  2(b)  we  attach  a mode  converter  module  to  the  sur- 
face of  the  waveguide.  As  wave  B strikes  the  module,  it  is 
converted  into  wave  C,  as  can  be  seen  in  the  decoupled  output. 
This  mode  converter  was  exposed  so  the  module  is  100*11.  ef- 
ficient in  converting  wave  B into  f. 

In  Fig.  2(c)  we  plage  a different  converter  module  on  the 


guide.  This  grating  has  index  modulation  such  that  half  the 
incident  energy  in  wave  B is  diffracted  into  wave  C.  This 
element  could  serve  as  a beam  splitter  for  use  as  the  first  ele- 
ment of  an  intrawaveguide  interferometer. 

It  should  be  noted  that  we  have  u,sed  a very  thick  optical 
waveguide  (6-mm  thick  Kodak  Microflat  glass).  Since  the 
guide  is  thick  and  the  surfaces  produce  a small  amount  of 
scattered  light,  we  can  observe  the  paths  of  individual  guided 
waves  (see  the  bright  spots  within  the  guide  in  Fig.  2).  For 
this  thick  guide  demonstration  we  also  do  not  have  to  worry 
about  waveguide  mrxles^  or  substrate  thickness  of  the  module. 
Our  modules  have  been  constructed  in  thick  dichromated 
gelatin  films  obtained  from  Kodak  649  F plates  ' Thus,  the 
emulsion  and  the  glass  substrate  (1  mm  thick)  of  this  film 
comprise  the  module.  Similar  devices  could  tie  built  for  thin 
waveguides. 

The  modular  optics  principle  can  be  extended  in  a number 
of  ways.  In  addition  to  constructing  the  modules  with  ex- 
ternal laser  beams  as  in  Fig.  1 , we  can  also  construct  modules 
by  exposing  a piece  of  sensitized  film  I which  has  been  tern 
porarily  attached  to  a guide)  at  the  intersection  point  of  the 
guided  outputs  from  two  other  modules.  We  also  need  not 
restrict  ourselves  to  plane  wave  devices,  as  mi«iular  lenses  and 
other  elements  can  easily  be  constructed.  .Active  elements 
can  be  included  to  build  up  modulators  and  switching 
a.ssemblies.  By  placing  a .series  of  different  modular  elemenus 
on  a waveguide,  entire  optical  circuits  can  be  assembled. 
Each  element  may  be  repositioned  and  adjusted  to  optimize 
the  performance  of  the  entire  waveguide  optical  circuit.  The 


Fig.  1,  A construci  ion  techniinie  lor  modular  t.pl  K ilt  elements.  4. 
H.  and  ( ’ are  incident  plane  waves;  /*,'  is  the  t hn  k flichromated  gi-lal  in 
emulsion  on  the  mialule  sulisirale  .S.  /*  is  an  ahsorplive  liacking  plate 
Xylene  is  used  as  an  index  matching  lluid  betwisui  the  components 
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elements  may  then  be  fixed  in  place.  As  it  is  presently  diffi- 
cult to  construct  an  entire  optical  circuit  via  an  integrated 
optics  approach,  the  modular  optics  building-block  concept 
appears  to  be  quite  useful. 

The  authors  wish  to  thank  Emmett  Leith  for  discussions 
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Fig.  2.  Modular  optical  elements:  (at  a coupler  module  is  used  to 
excite  guided  wave  H:  (h)  a mode  converter  module  converts  gtiided 
wave  H into  guided  wave  C,  (cl  a beam  splitter  module  converts  half 
the  energc  in  wave  K into  wave  C. 
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Directional  Coupling  Px'opariies  of  Interferorietric  VJavGguideG 

by 

G.  C.  Kung  and  W.  Y.  Vang 

Dn-3ctioaal  coupling  properties  of  sinusoidal  grating 
waveguides  are  investigated.  Our  conputer  analysis  which 
based  solely  on  the  diffraction  theory  of  thin  gratings 
predicl,s  results  that  are  in  good  agreement  with  those  calculated 
in  terms  of  coupled  mode  theory.  The  relationships  between 
directional  coupling  efficiency  and  some  parrsmeters  such  as 
grating  period,  input  pul.sQ  shape  and  guiding  channel  length 
are  disc'jsscd.  Our  analysis  predicts  that  we  can  substantially 
switch  the  input  energy  to  adjacent  channels  with  the  proper 
choice  of  parameters. 


Whan  this  work  was  done,  all  authors  were  with  the  University  of 
Ilichigan, Electrical  and  Conputer  Engineering  Department,  Ann 
Arbor,  Mi  UBlOh^  G.C.  Kung  is  currently  with  General  Dynamics, 
Pomona  Division,  .“omona , Ca  91766. 
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Introduction 


The  thin  grating  decomposition  (TGD)  method,  developed 

by  Alferness^  and  adapted  to  the  analysis  of  guided  wave 

2 3 

propogation  in  parallel-channel  planar  waveguides  ' (Leith 
et  al,  Rung  & Leith),  has  provided  a very  usual,  intuitive, 
and  simple  basis  for  examining  the  propagation  of  spatial 
pulses  of  arbitrary  shape  optical  waveguides  of  any  arbi- 
trary index  cross-sectional  profile.  Insights  are  provided 
in  situation  where  other  techniques  either  are  inapplicable 
or  extremely  difficult  to  apply. 

These  guides  can  alternatively  be  regarded  as  volume 
phase  gratings  of  various  index  perturbation;  the  wave  is 
guided  along  the  contours  of  maximum  refractive  index.  We 
consider  grating  of  sinusoidal  cross  section,  of  the  form 

t(x,z)  = exp  [jk  cos2TTx/d]  = 1 + jk  cos2iTx/d  (1) 

where  d is  the  grating  period. 

The  ability  to  confine  energy  within  each  guiding 
channel,  or  equivalently,  the  susceptibility  to  the  cou- 
pling of  light  energy  into  adjacent  channels  is  one  of  the 
most  important  characteristics  of  the  array.  We  use  the 
term  guiding  efficiency  (Hq)  to  mean  the  fraction  of  energy 
that  remains  in  the  input  channel  after  an  input  pulse  of 
energy  propagates  through  the  structure,  and  the  term  cou- 
pling efficiency  (rij^)  to  mean  the  fraction  of  energy  that 


26 


is  coupled  from  the  to  the  channel  over  the  inter- 

action distance. 

We  consider  the  propagation  of  parallel-channel  guides 
of  sinusoidal  profile,  and  we  consider  pulses  primarily  of 
the  form  of  a simple  cycle  of  a sinusoid.  We  report  various 
observation  which  show  the  power  of  the  technique  in  analogy 
a situation  not  easily  treated  other  methods,  and  therefore 
not  studied. 

Effect  of  Guide  Spacing 

We  examined  the  guiding  and  coupling  efficiencies  of 
guides  with  sinusoidal  index  variation  spacing.  The  input 
pulse  was  introduced  centered  at  a region  of  maximum  index 
and  has  the  form 

f(x)  = cos^  {—)  rect  (x/d)  (2) 

d 

where  d is  the  pulse  width.  This  pulse  thus  has  the  same 
shape  and  width  as  the  guide. 

The  results  are  shown  in  Figures  1-3,  we  note  as  ex- 
pected, that  the  guiding  efficiency  (n^)  in  the  input 
channel  decreases  monotonically  with  increasing  quadrature, 
as  plotted  in  Ref.  3. 

Figure  2 shows  the  same  informations  plotted  with 
guiding  efficiency  as  a function  of  index  modulation  for 
various  grating  frequencies.  In  all  cases,  the  input  pulse 
is  adjusted  so  as  to  be  of  the  same  width  as  the  guide. 


n , ( X I 0 ) 
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Guiding  efficiency  versus  index  nr^odulation 
for  sinusoidal  grating  waveguides  of 
various  spatial  frecuencies ; (a)  100  1/fjn; 

(b)  125  1/ntn;  (c)  143  1/ir.t;  (d)  1G7  1/rra; 

(e)  200  l/:'Mi.  •’’rating  thickness  = 0.54  rnm. 
Input  pulse  wjililis  are  equal  to  grating 
period. 


Also,  we  show  (Fig.  3)  the  coupling  efficiency  is  a 
function  of  n^. 

When  there  is  little  coupling  to  non-ad jacent  channels, 
the  increase  of  energy  confinement  in  the  input  channel  with 
increasing  n^^  corresponds  to  the  continuous  decrease  of 
energy  coupled  to  adjacent  channels  (curves  a and  b of  Fig. 
3) . However,  for  higher  frequency  gratings  such  as  curves 
c,  d and  e of  Fig.  3,  actually  increases  with  n^  in  the 
region  of  low  index  modulation,  this  occurs  because  a sig- 
nificant amount  of  energy  in  a grating  waveguide  of  high 
spatial  frequency  is  coupled  to  the  non-adjacent  channels 
when  the  index  modulation  is  small.  With  increasing  n^^ , the 
coupling  to  non-adjacent  channels  decreases.  Therefore 
results  in  less  leakage  from  the  input  channel. 

Effect  of  Input  Pulse  Shape 

The  Alferness  thin  grating  decomposition  analysis  is 
not  restricted  to  any  specific  profile  of  input  pulse.  We 
have  chosen,  for  simplicity,  a sinusoidal  pusle,  as  given  by 
Eq.  2.  Although  sinusoidal  grating  waveguides  cannot  pre- 
serve the  phase  information  (and  hence  the  shape)  of  the 
input  signal,  we  find  that  if  the  input  pulse  width  is  less 
than  or  equal  to  the  grating  period,  then,  with  increasing 
nj^,  the  output  pulse  width  converges  to  that  of  the  input 
pulse.  This  phenomenon  can  be  understood  by  consideration 
of  the  mode  structure  of  the  sinusoidal  grating  waveguides. 
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Coupling  efficiency  to  adjacent  channel 
vs.  index  modulation  for  sinusoidal  gratin 
waveguides  of  various  spatial  frequencies: 
(a)  100  1/mm;  (b)  125  1/mm;  (c)  143  1/mm; 
(d)  167  1/mjn;  (e)  200  1/mmi.  Grating  thick 
ness  = 0.54  mm.  Input  pulse  v/idths  are 
equal  to  grating  period. 


In  general,  when  a signal  is  launched  into  a waveguide, 
it  is  coupled  to  the  waveguide  modes,  some  which  are  guided 
modes  but  most  are  not.  The  transmission  of  the  input  pulse 
depends  entirely  on  the  guided  modes  which  are  excited  by 
the  input  field.  What  we  observe  at  the  output  is  a field 
distribution  which  is  a superposition  of  the  fields  in  the 
various  guided  modes;  in  most  instances  this  distribution 
bears  no  resemblance  to  the  shape  of  the  input  signal.  Of 
course,  the  input  coupling  is  most  efficient  when  the  input 
signal  is  matched  with  one  of  the  guided  mode  distributions; 
otherwise,  a part  of  the  energy  either  becomes  evanescent  or 
simply  propagates  through  the  structure  unguided. 

For  a sinusoidal  grating,  the  solution  of  the  wave 
equation  for  the  TE  modes  has  the  form 


Ej^(x,z)  = exp(-jB^z)  F(x)  (3) 

with  F(x)  satisfying 
d^F 

— j + (a  - 2q  cos2w)F  = 0 (4) 

dw'^ 

where 

W = TTfx  (5) 


Eq.  4 is  a Mathieu's  equation  whose  solution  may  be  ex- 
pressed in  the  form 

F (w)  = exp  (pw)j^Xcx,  By^exp  ( j2hv/)  (8) 
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where  p is  a constant.  We  consider  a specific  example  of  a 

_3 

grating  with  nQ  = 1.49  and  n^^  = .6  x 10  . From  Eg.  7 q is 
found  to  be  -0.893  when  Xq  is  chosen  to  be  0.6328  um.  The 
relation  between  a,  q,  and  p can  be  illustrated  in  the  form 
of  a stability  diagram  (Fig.  4) , following  the  standard 
treatment  of  Mathieu's  equation. ^ The  a-q  plane  is  divided 
into  various  stable  and  unstable  regions;  the  unstable 
solutions  correspond  to  complex  values  of  p,  which,  are  the 
guided  modes  of  the  structure.  The  operating  line,  q = 
-0.893  passes  through  several  unstable  regions.  From  the 
corresponding  values  of  a and  with  the  aid  of  Eq.  5,  we  can 
determine  the  z-directed  propagation  constant  B of  the 
guided  modes  whose  transverse  distributions  have  the  form  of 
Eq.  8. 

From  the  a-q  diagram,  it  is  evident  that  for  a fixed 
value  of  q,  there  are  many  continuous  sets  of  "a"  which 
correspond  to  many  continuous  sets  of  guided  modes.  How- 
ever, for  small  values  of  q,  the  ranges  of  the  guided  modes 
are  rather  narrow.  For  the  present  case  of  q = -0.893,  the 
values  of  "a"  which  give  rise  to  guided  modes  lie  in  the 
following  ranges:  a < -0.369,  0.018  < a < 1.802,  3.934  < 

a < 4.304,  9.039  < a < 9.061  and  so  on.  The  corresponding 
ranges  of  B/  obtained  from  Eq.  5 are:  B > 14.7959,  14.7944 

> B > 14.7885,  14.7814  > B >14.7801,  14,7643  > 8 > 14.7642 
and  so  on.  Thus,  even  though  the  guided  modes  in  theory 
form  many  continuous  sets,  the  ranges  of  B»  for  small 
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values  of  q,  are  so  narrow  that  the  modes  of  an  interfero- 
metrically  formed  waveguide  with  index  modulation  typically 
in  the  order  of  10  ^ can  be  considered  to  be  almost  discrete. 

This  would  explain  why  the  results  of  our  computer  analysis 
of  a sinusoidal  grating  waveguide  are  in  good  agreement  with 
those  obtained  from  the  conventional  analysis  of  single  mode 
coupled  waveguides  (i.e.,  those  with  discrete  boundaries). 

Pursuing  the  question  of  why  a sinusoidal  grating  wave- 
guide can  preserve  the  shape  of  an  input  pulse  of  the  form  f 

of  Eq.  1,  we  consider  a specific  example  by  arbitrarily 
choosing  a guided  mode  with  a value  of  a = 1.5;  resulting  in 
a p value  of  approximately  +0.3345  + j.  The  transverse  dis-  ■ 

tribution  of  this  mode  consequently  has  the  form 

F(x)  = exp  (-0. 33457Tfx^^^  Bj^exp  ( j2TTfhx) 

= exp  (-0 . 105x)  Q 2Bj^cos2iTfkx  (9) 

which  is  a set  of  slowly  damped  standing  waves.  For  an  in- 
put pulse  with  the  form  of  Eq.  9 and  equal  in  width  to  the 
grating  period,  i.e.,  d = 1/f,  we  immediately  see  that  the 
input  pulse  matches  closely  the  fundamental  space  harmonic 
of  the  guided  mode  except  for  the  slowly  damped  function, 
exp (-0 . 105x) . Consequently,  this  is  the  mode  which  is  prin- 
cipally excited;  what  we  observe  at  the  output  plane  is 
nothing  more  than  this  mode  distribution;  hence,  the  output 
pulse  converges  in  shape  to  the  input  pulse.  To  further 
illustrate  this  point,  we  propagate  a spatial  pulse  of  a 
quite  arbitrary  profile  through  a sinusoidal  grating  of 

Jj 


55 


thickness  0.54  irm.  The  distributions  of  the  input  and  out- 
put pulses  are  shown  in  Fig. 5.  From  these  results  and  also 
those  shown  in  Fig.  5 of  Ref.  3 for  an  input  pulse  extended 
over  three  channels,  it  is  quite  evident  that  the  output 
pulse  distribution  always  has  a nearly  co-sinusoidal  form 
within  the  guiding  channels,  regardless  of  what  form  the  in- 
put pulse  takes. 

Effect  of  Channel  Length 

The  continuous  exchange  of  energy  among  guiding 

channels  is  illustrated  in  Fig.  6.  An  input  pulse  of  3.01 

mm  width  is  placed  at  the  center  of  a guiding  channel  of  the 

same  width.  The  refractive  index  modulation  of  this  thick 

-4 

grating  is  arbitrarily  chosen  to  be  6 x 10  . As  the  pulse 

propagates  through  the  structure,  energy  of  the  pulse  is 
continuously  coupled  to  the  adjacent  channels.  The  results 
shown  in  Pig.  6 are  in  qualitative  agreement  with  that  ob- 
tained from  an  optical  discretely  bounded  single  mode  di- 

5 

rectional  coupler.  Using  the  quantity  guiding  efficiency 
and  the  same  parameters  as  before,  we  plot  the  fraction  of 
energy  in  the  0 and  +1  channels  against  propagation  depths. 
The  results  (Fig.  7)  show  energy  continuously  leaking  out  to 
the  adjacent  channels.  The  close  resemblance  between  inter- 
channelized  waveguides  and  directional  couplers  suggests  a 
quantitative  comparison  of  the  two;  although  there  are 
differences:  one  is  a single  mode  device  and  the  other  is 
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^ Propagation  of  an  input  pulse  of  arbitrary 
shape  in  a simn  oidal  grating  waveguide  of 
thickness  0.54  iiui,  l.iequeney  100  1/imn,  index 
modulation  .0  lO"^.  (a)  sliows  the  input 

pulse;  (b)  iihows  rhe  output  pulr;.'. 
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PROPAGATION  DEPTH  (mm) 


^ Fig.  7 Directional  coupling  of  grating  waveguices: 

, . (a)  shov;s  the  energy  that  remains  in  the 

input  channel  at  various  propagation  depths; 
(b)  shows  a function  Jo^  (.684s);  (c)  shows 

the  energy  directionally  coupled  to  t.he 
adjacent  channel;  (d)  shov;s  a fvu\ction 
Jl^  (.6842) 
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multi-mode  and  in  one  there  are  two  guides  in  paralles, 
whereas  in  the  other  there  are  many.  The  wave  amplitudes, 


aj^(z),  in  various  channels  of  the  single  mode  coupled  wave- 
guide when  the  energy  is  incident  on  the  k = 0 channel  is^ 

aj^(z)  = (-1)  (2Cz)  exp  (-j  6k^)  (10) 

where  J,  is  the  Bessel  function  of  the  first  kind,  C is  the 
' coupling  coefficient,  z is  the  propagation  depth  and  k is  an 

integer  denoting  the  channel  numbers.  The  energy  in  various 
channels  is 

Ej^(z)  = aj^^(z)  « (11) 

which  indicates  that  the  energy  in  the  0th  channel  varies  as 

the  square  of  a zero  order  Bessel  function.  We  superimpose 
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on  Fig.  7 the  functions  Jq  (.684z)  and  (.684z).  Note 
' that,  except  for  a small  region  of  z < 0.5  mm,  the  computed 

results  for  the  grating  waveguides  follow  very  closely  the 
two  Bessel  functions;  the  coupling  coefficient,  C,  of  the 
grating  waveguides  is  consequently  found  to  be  equal  to 
! 0.342. 

Optical  Directional  Coupler  Application 

i The  results  shown  in  Fig.  6,  predict  the  directional 

1 coupling  properties  of  a sinusoidally  stratified  structure; 

these  results  agree  qualitatively  to  those  calculated  in 
terms  of  coupled  mode  theory.  Some  additional  results  il- 
I ustrating  the  phase  relation  of  the  energy  between  the 


to 


coupled  guides  are  shown  in  Fig.  8.  By  placing  the  input 
pulse  of  0.01  Iran  width  in  front  of  the  entrance  surface  of 
the  guiding  structure  so  that  the  pulse  spreads  over  three 
channels  (Fig.  8a) , the  light  entering  into  the  side 
channels  suffers  a time  lag.  Our  analysis  predicts  (for 
properly  chosen  parameters)  constructive  interference  of  the 
light  entering  directly  into  the  side  channels  and  that 
directionally  coupled  into  the  side  channels  from  the 
central  one,  because  the  directionally  coupled  light  also 
suffers  a time  lag.^  On  the  other  hand,  by  placing  the  in- 
put pulse  sufficiently  behind  the  entrance  surface  (Fig.  8b), 
the  interference  in  the  side  channels  is  destructive  because 
now  the  direct  lights  entering  the  side  channels  have  a time 
lead.  Somekh  has  experimentally  demonstrated  such  phenom- 
enon. It  is  noteworthy  that  our  analysis,  which  is  based  on 
nothing  more  than  a continuous  iteration  of  the  quite  simple 
thin  grating  diffraction  theory,  should  so  readily  predict 
these  rather  subtle  effects. 

Our  analysis  indicates  that  we  can  substantially  shift 
the  input  energy  to  adjacent  channels  of  multichannel  wave- 
guides by  using  the  proper  choice  of  parameters.  We  can 
envision  various  applications,  such  as  incorporation  of  this 
feature  into  an  acousto-optical  divice  using  the  modulated 
medium  as  the  guiding  channels,  which  is  possible  provided 
sufficient  index  modulation  is  achieved. 

Since  the  index  modulation  can  be  controlled,  it  would 
be  possible  to  efficiently  switch  the  input  energy  between 
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Phase?  relation  betv/een  directionally 
coupled  modes  in  a sinusoidal  grating 
waveguide.  To  the  left  are  sketches 
illustratiinr  varieiis  input  conditions. 
The  pcsitioiui  ej;  the  input  pulse  ,'ire: 

(a)  0.21  mm  in  J i;ont  of  the  grating; 

(b)  0.27  mm  ’.behind  the  front  surf. ice  of 
the  grating;  (c)  on  the  front  suriace. 

To  the  right  arc  the  correspendiner  in- 
tensity distributions  at  the  output 
plane.  Thickness  of  grating  is  2.16  mm, 
grating  frequency  is  100  1/mm  and  index 
modulat4.on  is  0.6  x lO”^. 


r ^ 

^ the  0th  and  1st  channels,  thus  forming  a simple,  yet  versa- 

[ tile  multi-mode  switching  device. 
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Abstract 

A waveguide  interferometer  is  assembled  using  holographi- 
cally formed  optical  elements  that  are  constructed  as  individual 
components  and  can  be  selectively  attached  to  an  optical  wave- 
guide. The  interferometer  is  very  rugged  and  easy  to  align. 


'Work  performed  at  The  University  of  Michigan. 
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1.  Modular  Optical  Elements 

We  describe  modular  optical  elements  as  miniature  optical 

components  that  are  individually  constructed  and  can  be  selectively 
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attached  to  the  surface  of  an  optical  waveguide.  ’ ’ These 
components  can  be  Individually  adjusted  to  optimize  system  perfor- 
mance. In  this  paper,  we  describe  the  construction  and  use  of 
modular  optical  elements  to  form  a waveguide  interferometer. 

Our  optical  modules  consist  of  holographic  optical  elements 
constructed  in  dichromated  gelatin  films  that  are  coated  on  glass 
substrates.  Each  element  serves  as  both  a grating  waveguide- 
coupler  and  a beam  splitter,  as  shown  in  Fig.  1.  The  modular 
gratings  are  constructed  by  incoherently  superimposing  two  thick 
gratings  within  a dichromated  gelatin  film  such  that  the  two  grat- 
ings have  a common  Bragg  angle,  A.^  Readout  at  the  angle  A will 
produce  two  diffracted  waves  each  of  which  contains  half  of  the 
incident  energy.  The  two  diffracted  waves  are  at  angles  such 
that  they  will  both  be  trapped  within  the  guide  G. 

2.  The  Interferometer 

Having  constructed  our  grating  modules,  we  assemble  our  in- 
line interf cjrometer  by  placing  two  of  these  modular  gratings  on  the 
surface  of  uhe  thick  glass  optical  waveguide  as  shown  in  Fig.  2. 

Xylene  is  used  as  an  index  matching  liquid  between  the  modules  and 
the  guide.  The  first  element  (on  the  left)  is  used  as  the  waveguide 

coupler  and  beam  splitter  while  the  second  element  acts  as  a beam-  i 


combiner  and  decoupler. 


The  vavepuide  interferometer.  Proceeding 
to  right,  one  can  see  the  coupler  module 
’’bject,  tlie  decoupler  module  and  the  out- 
•.•;hich  strikes  the  white  card. 


L 


The  Interferometer  is  used  by  illuminating  the  coupler  module 

with  a coherent  beam  of  light  from  the  upper  left  in  Fig.  2.  Since 

the  incident  wave  (A  in  Fig.  1)  has  a finite  cross  sectional  area 
2 

(=  1 cm  ) , the  waveguide  is  thick  (=  6 mm)  and  the  waves  B and  C 
propagate  at  greatly  different  angles  (=  20°  difference),  the  areas 
at  which  the  guided  waves  strike  the  surface  of  the  guide  are 
physically  separated.  Thus,  we  can  insert  a sample  into  the  inter- 
ferometer by  laying  it  on  the  surface  of  the  waveguide  at  a position 
where  only  one  cf  the  guided  waves  (e.g.  wave  C)  strikes  the  surface. 
An  index  matching  liquid  (in  this  cas^  xylene)  is  placed  between  the 
sample  and  the  guide.  Thus,  guided  wave  C leaves  the  guide,  passes 
through  the  sample,  and  re-enters  the  guide  as  a modulated  wave. 

A short  distance  later,  where  waves  B and  C again  overlap,  the  beam 
combiner  has  been  placed  on  the  surface  of  the  guide.  This  modular 
element  in  identical  in  construction  to  the  beam  splitter  element  and 
is  merely  rotated  by  180°  (about  a vertical  axis  in  Fig.  1)  before 
being  placed  on  the  guide.  The  two  guided  waves  simultaneously  enter 
the  final  grating,  are  combined,  and  are  diffracted  out  of  the  guide. ^ 
The  combined  waves  leave  the  guide  at  approximately  45°  as  can  be 
seen  in  Fig.  2 where  they  strike  a white  card  after  leaving  the  inter- 
ferometer . 

3.  Experimental  Results 

The  output  of  the  interferometer  is  photographer  by  placing  a 
camera  in  place  of  the  white  card  in  Fig.  2.  With  nt  sample  in  the 
interferometer,  we  see  straight  fringes  present  on  the  output  (Fig.  3). 
The  fringe  frequency  in  Fig.  3 can  easily  be  varied  by  small  rotations 

1 Q 


Fig.  3.  The  interferometer  output  showing  straight 
fringes  when  no  test  object  is  in  the  interferometer. 


Fig.  A,  Interferometer  output  when  a piece  of  thin 
glass  with  a deep  scratch  is  used  as  a test  object. 


!0 
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of  the  final  grating  (about  a vertical  axis  in  Fig.  1).  Because 
the  angles  between  guided  waves  B and  C are  fixed  during  grating 
construction  and  cannot  get  out  of  alignment,  the  only  adjust- 

4 

ment  that  ever  needs  to  be  done  is  a slight  rotation  of  the 
gratings. 

In  Fig.  4,  we  see  the  output  when  an  object  has  been  insert- 
ed into  the  Interferometer.  The  test  object  is  a small  Tiece  of 
glass  (microscope  slide  cover-glass)  in  which  we  have  put  a deep 
scratch  with  a glass  cutter.  The  scratch  extends  half  way  across 
our  test  object  and  shows  up  as  dislocated  fringes  on  the  left 
side  of  the  output. 

4.  Conclusions 

The  modular  interferometer  we  have  demonstrated  is  easy  to 
adjust  and  is  very  rugged.  Rapid  changes  in  air  temperature  or 
air  currents  have  little  effect  on  the  device. 

The  modular  coupler  elements  can  easily  and  quickly  be 
placed  on  any  piece  of  glass  which  we  would  like  to  interfero- 
metrically  test  for  flatness  or  other  optical  properties. 
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Abstract 

A method  of  analyzing  optical  propagation  in  thick  holographic  reflection 
grating  by  decomposition  of  the  thick  grating  into  thin  slabs  is  investigated. 
The  method  is  applicable  to  study  higher  order  Bragg  diffraction  process  in 
thick  reflection  grating.  The  method  is  shown  analytically  to  be  equivalent 
to  the  coupled-wave  solution  of  Maxwell's  equations  for  a thick  sinusoidal 
reflection  grating. 


1.  Introduction 


A physically  intuitive  method  for  analyzing  wave  propagation  in  thick 
holographic  transmission  grating  was  developed  by  R.  Alferness^  in  which  the 
thick  gratings  are  mathematically  decomposed  into  a series  of  thin  slabs,  with 
each  acting  as  a thin  grating.  For  a readout  plane  wave  of  arbitrary  spatial 
frequency,  the  total  amplitude  of  the  zero  and  dlssracted  orders  is  derived 
from  computing  the  succesive  effect  of  each  thin  grating  is  determined  solely 
by  thin  grating  theory.  The  thin  grating  decomposition  method  (TGD)  for 
transmission  grating  has  been  widely  used  to  study  higher  order  diffraction 
and  waveguide  phenomenon.  However,  the  coupling  matrix  does  not  apply  to  the 
case  of  reflection  gratings.  Recently,  the  use  of  reflection  gratings  has 
been  considered  for  Integrated,  in  such  application  as  wavelength  fibers  and 
distributed  feedback  (DFB)  and  distributed  Bragg  reflector  (DBR)  lasers.  In 
this  paper,  the  TGD  coupling  matrix  method  is  extended  to  the  reflection 
grating  case. 

Our  method  for  analyzing  optical  propagation  in  volume  reflection  gratii.^ 

offers,  because  of  its  conceptual  simplicity  and  ease  of  implementation, 
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advantages  over  the  theory  given  by  Burckhart  and  Kasper  . Like  the  matrix 
theory  the  method  is  numerical;  it  requires  only  matrix  multiplication  and 
it  proceeds  to  decompose  the  thick  reflection  grating  into  a series  of  thin 
slabs,  each  of  which  acts  simply  as  a thin  grating. 


2.  Description  of  the  Method.  We  describe  the  method  with  a single  re- 
flection grating  structure  produced  holographically,  as  shown  in  Fig.  1. 

The  film  emulsion  of  thickness  of  D is  perpendicular  to  the  z-axis  and  is 
assumed  to  be  of  infinite  extent  in  the  x and  y directions.  Plane  wave  R 
of  amplitude  a^  incident  at  angle  0^^  (angle  measured  inside  the  emulsion) 
with  respect  to  the  z-axis,  and  plane  wave  S of  amplitude  at  6^  interfere 
to  produce  a grating  throughout  the  volume  of  the  emulsion.  The  exposure 
E = IT  (I  is  the  total  intensity  and  T is  the  exposure  time)  can  be  written  as 

E(x,z)  = ^ cos  2iTf(z  - z^)]  , 

where  f = n^(cos  9^^  + cos  6^)/A  is  defined  as  the  spatial  frequency  in  the 
z-direction,  X is  the  free  space  wavelength,  z = x tan  6,6=  (9  - 6„)/2 

X 1 J 

2 2 '^2 

is  the  slant  of  the  fringes,  = a^  + ^2  ’ ^ ~ 2a^a2/(a“  + ^2^  ‘ 

any  plane  x = constant  with  Az  as  the  thickness  of  the  thin  slab,  the 

fringes  are  sinusoidally  modulated  with  frequency  f.  However,  because  of 

fringe  slant  (when  the  construction  geometry  is  not  symmetrical  about  the 

x-axis)  the  fringes  at  z are  shifted  by  z^  relative  to  those  at  x = 0, 

Introduce  a relative  phase  shift  of  -2iTfz  . 

x 


which 


We  examine  the  properties  of  optical  propagation  in  the  thick 

reflection  grating  of  Fig.  1 by  decomposing  it  into  a series  of  thin 

slcdas  of  thickness  Az  which  must  be  sufficiently  small.  The  proper 

upper  limit  on  Az  is  that  for  a fixed  readout  wavelength  and  grating 
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frequency,  Kogelnik's  Q factor  (Q  = 2'iTXAzf'^/n^ ) for  the  thin  grating 
be  less  than  one.  Each  of  the  thin  gratings,  therefore,  exhibits  no 
Bragg  effects  and  its  angular  distribution  of  diffracted  light  is  cal- 
culated based  on  the  thin  grating  theory.  We  assume  a pure  phase 
reflection  hologram  in  which  the  diffraction  process  takes  place  dis- 
cretely at  the  plane  z = nAz,  n = 1,2,...  . We  also  assume  a linear 
relationslup  between  the  index  of  refraction  and  the  exposure 
E:  n(x)  = ME.  Therefore,  the  amplitude  transmittance  as  a result  from 
the  exposure  given  by  Eq.  1 for  a plane  wave  at  0^  Is  expanded  according 


to  the  method  by  Alferness  : 


where  J (b)  is  the  order  Bessel  function  b = (2AAzn,/^  cos  9 ) 
q 1 1 

and  z^=x  tan  6.  The  n^^  thin  grating  can  be  represented  by  an  infinite 
superposition  of  amplitude  gratings  of  spatial  frequency  qf  with  diffrac- 
tion amplitudes  equal  to  the  qth  order  Bessel  function.  The  terms 
with  minus  q portion  in  the  summation  sign  correspond  to  the 
reflected  diffracted  waves  and  the  terms  with  positive  q portion  corres- 
pond to  the  transmitted  diffracted  wave.  The  exponential  outside  the 

summation  sign  in  Eq.  3 represents  an  unimportant  common  factor  and 

th 

can  thus  be  neglected.  The  q frequency  component  is  shifted  in 
phase  by  -2qf  x tan  9 because  of  the  fringe  slant.  This  phase  deter- 
mines the  Bragg  angle  in  the  slanted  reflection  grating. 

To  find  the  effect  of  the  thin  slab  upon  an  incident,  a forward 
and  a backward  plane  wave  of  arbitrary  spatial  frequency,  we  extend 
the  thin  grating  theory  to  include  the  propagation  factors  of  both 
forward  and  backward  plane  waves  of  arbitrary  spatial  frequency  in  the 
thin  slab.  First,  propagation  ^-'Ctors  of  the  incident  waves  over  the 
interval  can  be  described  by  the  free  space  transfer  function.  Then, 
the  separation  of  incident  forward  or  backward  waves  into  the  various 
diffracted  waves  at  2 = nAz  is  described  by  the  amplitude  transmittance 
function  in  Eq.  3.  Therefore,  for  an  input  forward  plane  wave  of 
x-spatial  frequency  f.,  amplitude  A.^  and  phase  9^  at  z = (n  - l)Az, 
the  field  at  z = nAz  (i.e.,  out  of  the  n^*^  grating)  is 

+00 


Aj{x,nAz)  = exp  (j(})^)  ^ ( j ) J (b)  exp  (-  j27rqfz^) 


q=-o 


exp  (jf^x)  exp 


2irn 


0 r 


1 - Uf 


./n  ) Az 
10 


exp  (j2Tiqfz)  exp  ' 


(4) 


t 


where  f.  = f.  - 27iqf  tan  6 and 
iq  1 


f . = /I  - [ (A/n 

iZ  0 iq 


Similar  results  can  be  obtained  in  an  input  backward  plane  wave  except 
for  a minus  sign  in  the  directional  cosine.  Therefore,  each  plane  wave 
incident  upon  the  n^*^  slab  generates  a countable  number  of  plane  waves 
with  z-spatial  frequency  + qf)  • The  amplitude  of  the  q^^  wave  is 

the  product  of  the  input  amplitude  and  the  strength  of  the  q^^  grating 
frequency  component. 

The  phase  of  each  wave  is  given  by  three  factors.  The  term 

exp  (-j2iTqfz  ) in  Eq.  A.  represents  a phase  shift  resulting  from  the 
n 

parallel  shift  of  the  n*^  grating  relative  to  the  (n  - l)*"^  (a  result 
of  the  fringe  slant) . The  term 


-f^ 


2 h/2 


exp  . 1 


gives  the  phase  shift  due  to  propagation  through  the  distance  Az.  In 
addition,  because  a phase  material  is  used,  there  is  the  phase  factor 
(j)*^.  Then  each  plane  wave  output  from  the  n^^  slab  is  an  input  to 
the  (n  + 1)  slab. 

The  thick  grating  is  illuminated  at  z = 0 with  a plane  wave  of 
z-spatial  frequency  f^^and  an  amplitude  . We  represent  the  amplitudes 
of  each  of  the  possible  plane  waves  that  exit  from  the  n^^  grating  by 


, i integer  = 0 , 


) 


where  is  the  amplitude  of  the  forward  wave  of  z-spatial 

frequency  while  is  the  total  amplitude  of  the  backward 

wave  of  spatial  frequency  [f^  - (i  + l)f].  Because  Eq.  4 is  valid 
for  an  input  plane  wave  of  arbitrary  spatial  frequency,  we  can  relate 
the  total  amplitude  of  each  of  the  plane  waves  at  z = n^z  to  its  value 
at  z = (n  - l)Az  by 


HA  , 
n n- 1 


where  the  matrix  of  coupling  coefficients  H,  . determined  from  Eq . 2.19 

n k]  ^ 

describe  the  eimplitude  and  phase  shift  effected  by  the  n^^  thin  slab 
in  diffracting  the  plane  wave  into  the  k^^  plane  wave.  The  coupling 
matrix  for  the  pure  phase  case  appears  in  Appendix  A. 

If  we  use  the  initial  condition  at  z = 0 to  write 


the  amplitude  at  z = nAz  will  be 


A = H HA 
n n 1 o 


The  amplitude  of  each  of  the  plane  waves  at  the  exit  plane  z = D is 
given  by  Eq.  8 for  n = N = D/Az. 

Results  and  Their  Equivalence  to  Coupled  Wave 


Approach.  For  a grating  with  sufficiently  high  Q factor  it  is  assumed 
that  only  the  zero  and  one  diffracted  order  are  significant  for  readout 


about  the  first  order  Bragg  angle.  Because  only  two  orders  are 


important,  the  coupling  matrix  for  the  n 


thin  grating  in  Eq.  4 


can  be  truncated  to  a 2 x 2 matrix: 


e J (b  ) 
o o 


• ^ 36  (n) 

36  Jj(b^)e 


■je  (b)e-^«<"^ 

- 1 o 


e J (b  ) 
o 1 


where  b^  = 2TiAzn^/X  cos  0^, 

= 2ttAz  cos 
6 = 2irnfx  tan  6 , 

cos  6^  = [1  - (X/n^)2f2]^/2  ^ ^ ^ 0,1,...  , (10) 

= the  X -spatial  frequency  of  the  readout  wave, 

\ = the  free  space  wavelength  of  the  readout  wave, 

9 = the  slope  of  the  fringes  with  respect  to  the  x axis,  and 
n^  = the  index  modulation  of  the  grating. 

To  simplify  we  assume  that  the  fringe  slope  angle  0 is  zero  (i.e., 
unslanted  fringes) . For  typical  values  of  index  modulation  and  az 
arbitrarily  small,  we  can  approximate  the  Bessel  functions  by  their 
first  order  approximations. 


J (x)  = 1 J (x)  = x/2  for  X <<  1 

o 1 


The  coupling  matrix  for  each  thin  slab  can  therefore  be  written 


P 5 Az  cos  9 

exp  IX  o 


Jr 

cos  6 


-j  e^^*^ 

3 COS  9 J 


where  \ = 


rn  Az/A  a:i:i 
1 


2$  = (2TTnQAz/A)  (-  cos  0^  + cos  9^)/ 


(12) 


in  which  J ^^(x)  (x)  . All  angles  are  internal. 


2$  = 2TTn„Az/X  or  2<t  = -2nn  Az(AX)  cos  0 /A 

0 0 o 


(13) 


th 


If  the  amplitudes  of  the  zero  and  first  order  at  the  n thin  slab  are 
expressed  by  the  elements  of  a column  vector 


A(n)  = 


(14) 


then  the  value  of  A at  the  exit  plane  of  the  emulsion,  z = D is 


N. 


A(N)  = H A(0)  , 


(15) 


where  N = D/Az  and  initial  conditions 

at  z = D , R(0)  = R , S(0)  = 0 A(0) 


(16) 


and 


at  2 = 0 


R(N)  = 1 


S(N)  = S 


N 


A(N)  = 


(16) 


can  be  obtained  by  performing  the  matrix  multiplication  of  Eq.15  in 
the  limit  N -*■  «>,  Az  -*■  0 such  that  NAz  = D. 

As  shown  by  the  coupling  matrix  of  Eq.  11  2$  represents  a 
relative  phase  shift  between  the  two  forward  and  nackward  waves  that 
were  produced  by  a thin  slaib  when  the  reconstruction  wave  length  or 
angle  differ  from  their  construction  values.  It  is  this  phase  shift 
that  causes  the  decrease  in  diffraction  efficiency  for  readout  waves 
incident  off  the  Bragg  anqle;  4>  corresponds  to  the  dephasing  pareimeter 
of  the  couplcd-wave  theory. 


The  matrix  multiplication  of  Eq.  15  is  greatly  simplified 


when  H is  represented  in  a basis  such  that  it  is  diagonal.  The  diagonal 
elements  are  just  the  eigenvalues  of  H and  the  basis  vectors  are  the 
associated  eigenvectors.  The  eigenvalue  of  H can  be  easily  shown  to  be 


, = exp  (j$)  exp  ± 

^ 9 ^ 


where  we  have  used  the  limit  Az  0 to  set  sin^  6=6^.  Also,  C = -C 

o 1 

= cos  0 . The  associated  eigenvector  E , E are  then 
o 12 


1 ^2  ^1 


where  a.  = j (C, /y ) (1  - <» . ) e 
1 1 1. 

From  Eq.  11  and  the  initial  condition  on  A (Eq.  16),  we  find 


the  initial  condition  in  the  E basis;  E(0) 


diagonal  representation. 


. Therefore,  in  the 


E(N)  = 


R 


Rx 


Rx  ij). 


We  then  use  the  inverse  of  the  transformation  of  Eq. 18  to  transform 
back  to  tlie  original  representation  to  find  the  amplitude  of  the 
diffracted  order: 


Aj  (N) 


R , N ,N. 

Ta7— ,T  - V 


A (N) 
o 


Thus, 


-t(h 

A (N)  = Rxe-'  X 

o 


j*  sinh 


N2y2 


$2fj2 


1/2  r .2  ^1/2 

I $2 


cosh 


f 

MSS 


^2n2 


' 1/2'j 

1 

X 

J 

\ 

Jc:__  <^2 


c c 

IS  1 


■1/2 


Aj  (N)  = Rxj  X 


sinh 


nS2  22!^^^ 

Ji_. 


(22) 


C C 
0 1 


From  Alferness  ^ ,■  we  have  N^y^/C^C^  -*■  and  N^4’  -*■  C in  Kogelnik's 
notation,  therefore 


1/2 


• j* 

-je-’ 


A^  (N) 


jC/v  ■*-(!-  coth  (v^  - 


(23) 


which  is  identical  to  Kogelnik's  formula.  The  diffraction  efficiency  is 
C. 

DE(N)  = 1a  (N)]2  = l/[l  -t-  (1  - C^/v2)/sinh2  (v^  - . 

‘-0 

(24) 

We  have  thus  shown  analytically  the  equivalence  between  the  method  of 
reflecting  thin  grating  decomposition  and  Kogelnik's  coupled-wave  solu- 
tion for  wave  propagation  in  a thick  sinusoidal  reflecting  grating  for 
which  the  assumption  of  the  latter  theory  are  valid.  It  is  easy  to 
write  a computer  program  for  solving  any  nonuniformly  longitudinal 


D ( 


index  modulation  cases. 
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ABSTRACT 


We  present  an  approach  to  the  analysis  of  a holo- 
graphically produced  grating  coupler,  in  which  we  sepa- 
rate the  coupling  oroblern  into  independent  slab  wave- 
guide and  diffraction  grating  problems.  Considering  the 
guiding  structure  as  a slab  waveguide,  we  obtain  the 
paramenters  of  the  propagating  modes  of  the  coupler 
through  conventional  analysis.  The  diffraction  ef- 
ficiency of  the  grating  orders  are  then  calculated  by 
using  the  thin  grating  decomposition  method  (TGD) . 
Utilizing  the  propagating  mode  angles  of  the  waveguide 
and  the  calculated  diffraction  efficiencies  of  the 
grating  , we  are  able  to  calculate  the  radiation  loss 
coefficient  and  therefore  the  coupling  efficiency  of 
the  holoqraphic  waveguide  coupler. 
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The  basic  component  of  an  integrated  optical 
circuit  is  the  thin-film  dielectric  waveguide  upon 
which  a waveguide  coupler  is  built.  The  waveguide 
coupler  can  be  used  to  convert  an  incident  light 
beam  into  a mode,  or  modes,  of  the  waveguide. 
Couplers  which  are  integral  components  of  the  wave- 


guide are  either  corrugated  surface  gratings  or  holo- 
graphic gratings  above  the  waveguide.  We  suggest  a 
third  waveguide  coupler  configuration  in  which  the 
grating  is  fabricated  inside  the  waveguide  structure. 

We  initially  adapted  the  first-order  perturbation 

analysis  of  Ogawa^  to  this  structure  but  obtained  an 

excessively  high  radiation  loss  coefficient.  Peng  et. 

2 

al  used  the  rigorous  electromagnetic  boundary  problem 
to  analyze  the  holographic  grating  coupler  above  the 
waveguide,  but  this  method  seems  inapplicable  to  the 
case  where  the  grating  is  inside  the  waveguide.  We 
therefore  develope  a different  approach,  where  we 
separate  the  analysis  into  a diffraction  grating  analysis 
and  a slab  waveguide  analysis.  We  obtain  guided  and 
radiation  modes  as  solutions  of  the  unperturbed  wave- 
guide structure  using  a conventional  slab  waveguide 


^ 


analysis.  We  then  use  the  thin  grating  decomposition 
method  of  Alferness^  to  analyze  the  grating  diffrac- 
tion process  which  produces  an  exchange  of  energy 
among  the  guided  and  radiation  modes  as  the  guided 
wave  propagates  through  the  thickness  of  the  grating 
waveguide.  Combining  the  slab  waveguide  and  dif- 
fraction grating  results  we  then  determine  the  radi- 
ation loss  coefficient  by  accounting  for  the  energy 
lost  through  radiation  modes.  From  such  a description 
of  the  waveguide  coupler  we  then  determine  the  cor- 
responding coupling  efficiency. 
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Problem  Formulation 


The  coupler  problem  to  be  solved  is  described 
with  reference  to  Fig.l.  A plane  wave  (not  shown) 
impinges  on  the  grating  from  its  air  side,  at  an  angle 
which  satisfies  the  second  Bragg  condition  for  the 
grating  in  the  waveguide.  The  grating  is  constructed 
so  that  the  second  diffracted  order  is  a propagating 
mode  for  the  planar  waveguide  structure  containing  the 
grating.  As  this  diffracted  wave  propagates,  a portion 
of  its  energy  is  diffracted  into  the  first  and  second 
orders  thus  producing  a coupler  loss  mechanism. 

For  purposes  of  analysis  we  find  it  convenient  to 
reverse  the  process  using  the  concept  of  optical 
reciprocity.  We  assume  that  a guided,  propagating  wave 
(represented  in  Fig.l  by  Poynting  vector  P^^M  is 
incident  on  the  grating  portion  of  the  waveguide.  As 
it  travels  it  produces  diffracted  waves  P^  and  P^ 

By  calculating  the  magnitude  of  these  diffracted  waves, 
we  can  obtain  the  amount  of  incident  light  remaining 
in  p^*^^  and  lost  into  P^  and  P^  We  can  there- 

fore obtain  the  radiation  loss  coefficient  and  the 
coupling  efficiency  which  characterize  the 
waveguide  coupler. 


^1 


grating 


FIG.  I 

Schematic  diagram  of  grating  waveguide  configuration  with 

guided  mode  (zeroth  'rder)  Poynting  vector  and  radiation 

mode  (first  and  second  order)  Poynting  vectors  (P^  ^ 

i 

I 


From  a geometrical  optics  viewpoint  we  solve  the 

planar  waveguide  problem  for  the  mode  reflection  angle 

(0  ) , the  effective  guide  width  (W  ) and  the  Poynting 
m m 

vector  (P^*^^)and  alsodefine  a mode  bouncing  rate  B . 

m 

We  assume  that  the  guided  wave  with  amplitude  A is 
reflected  back  and  forth  between  the  waveguide  bounda- 
ries at  angles  and  -0  . The  effective  guide  width 
mm  ^ 

is  obtained  from  the  expression  W ~D+p^+p^, 

where  p^  and  p^  are  the  exponential  decay  rates  for 
a S 

mode  m in  the  air  and  substrate  regions  with  refrac- 
tive indices  n and  n respectively.  The  bounce  rate 

a S 

for  mode  m in  terms  of  0 and  W becomes, 

m m 

= [2W  tan0  • (1) 

m ‘ m m 

In  Fig.  1 there  are  two  diffracted  wave  orders  with 
Poynting  vectors  P^  and  P^  Using  these  Poynting 

vectors,  a diffraction  efficiency  can  be  defined  as 
the  fractional  amount  of  incident  power  diffracted 
into  the  various  orders.  Defining  these  ratios  of 
power  flow  in  the  x direction  in  terms  of  amplitude 
ratios, we  obtain 


COS0^^> 


COS0 

m 


where 


^(q) 


A. 

1 

0 

m 

^(q) 

e 


= diffraction  efficiency  of  order 

= amplitude  of  order  wave  after  propagation 

through  grating 

= amplitude  of  incident  wave  on  grating 


= angle  of  propagating  mode 


(q)  _ 


th 


= external  angle  of  diffracted  q order,  where 


n sin0  = n sin0 
g e e 


th 


internal  angle  of  diffracted  q order 
index  of  external  medium  of  q^*^  order 
bulk  index  of  grating  • 


It  should  be  noted  that  the  diffraction  efficiency  of 
a particular  mode  depends  upon  whether  it  remains  in 
the  waveguide  or  escapes  to  the  substrate  or  air. 

This  dependence  is  accounted  for  by  setting  n^  equal 
to  the  index  of  the  material  in  to  which  the  q^^  order 
escapes  or  is  trapped,  0^  being  the  external  angle  in 
this  medium.  By  using  the  TGD  method  we  can  determine 
the  amplitudes  of  the  different  Bragg  diffraction 
output  orders  for  a given  amplitude  input  (A^). 

The  analysis  of  conventional  bulk  phase  gratings 
is  applicable  to  the  thin  film  grating  coupler.  The 
basic  difference  between  these  two  grating  types  is 


'4 


that  a bulk  grating  diffracts  one  free  space  propa- 
gation beam  into  another,  whereas  the  Bragg  thin  film 
grating  coupler  transforms  a free  space  propagation 
mode  into  a leaky  guided  mode  of  the  waveguide.  The 
leaky  mode  m is  coupled  into  transmitted  radi- 
ation modes  via  diffraction  by  the  same  grating  orders 
which  coupled  the  light  into  that  leaky  auided 
mode . 

The  power  lost  over  the  length  of  the  grating 
coupler  is  characterized  by  the  power  loss  coefficient, 
or  equivalently  by  the  amplitude  loss  coefficient 
a . Because  this  loss  is  caused  by  diffraction  of 
the  guided  wave  by  the  grating  ,the  loss  coefficients 
become  directly  proportional  to  the  diffraction  ef- 
ficiency of  each  diffracted  order.  We  proceed  to 
obtain  the  characteristic  amplitude  loss  coefficient 
in  terms  of  the  diffraction  efficiencies  by  con- 
sidering a single  section  of  the  grating  parallel  to 
the  x-axis.  We  then  extend  the  results  over  the 
entire  length  of  the  grating  to  obtain  a radiation 
(amplitude)  loss  coefficient  and  the  coupling  ef- 
ficiency of  the  grating  coupler. 
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Loss  Mechanisms  and  the  Coupling  Coefficient 

In  the  following  approach  we  use  the  ray  optics 

method  to  obtain  an  expression  for  the  amplitude  of 
the  leaky  guided  mode  and  the  radiation  loss  coef- 
ficient. Referring  to  Fig. 2 we  divide  the  total 
grating  length  into  sections  (2W^tan0^)  over  which 
one  bounce  of  the  leaky  guided  mode  occurs.  We  note 
that  this  is  the  distance  in  the  z direction  separating 
equi-phase  points  along  the  ray  path.  We  qualify  the 
loss  over  this  section  as  the  decrease  in  incident 
amplitude  (A^)  diffracted  out  of  the  leaky  guided 
mode  as  it  satisfies  the  Bragg  diffraction  condition 
for  the  grating  over  one  half  of  the  section  length. 
Therefore  the  guided  wave  amplitude  after  propaga- 
tion through  one  section  is, 

A = 11-'.'^’-  (3) 

We  now  consider  the  field  at  some  plane  z = UB  ^ . 

m 


Allowing  the  guided  wave  to  propagate  through  £ 
sections  of  the  grating  coupler  and  normalizing  A^ 
to  one,  we  obtain  for  the  leaky  guided  mode  amplitude, 


FIG.  2 

One  unit  of  the  sectionalization  of  the  grating  structure, 

11/  9 ' / 

illustrating  the  decoupled  energv,  A^(n,)  ^ ^ and  A (ri,)  '^^with  the 

o a o d 

remaining  coupled  energy.  Loss  occurs  only  over  one  half  of 
the  section  where  the  guided  mode  satisfies  the  Bragg  condition. 


where  we  have  obtained  an  expression  for  the  leaky 
guided  mode  amplitude  along  the  length  of  the  coupler. 


We  define  the  total  radiation  loss  coefficient  (a  ) 
to  be , 


a = 
a 


4W  tan9 
m m 


The  leaky  mode  amplitude  therefore  decays  as 
exp  {-a  z)  for  z > 0.  Since  the  energy  lost  in  the 
leaky  guided  mode  must  be  gained  by  the  radiation 
mode,  the  radiation  mode  in  the  near  field  outside 
the  waveguide  must  have  an  exponential  shape  comple- 
mentary to  the  guided  leaky  mode.  The  amplitude  of 
the  radiation  mode  can  be  calculated  from  the  fact 
that  for  a semi-infinitely  long  grating  the  energy 
of  an  incident  guided  mode  at  z = 0 must  be  equal  to 
the  integrated  energy  in  the  radiation  mode.  The 
energy  becomes, 

1 = /"e^  (z)dx  (8) 

where  we  set  the  integral  equal  to  the  initial  normal- 
ized energy  of  one  in  the  guided  mode.  The  electric 
field  of  the  radiation  is  then  given  by, 


E (z)  = E exp{-a  z) 
r o a 


(9) 


Substituting  Eq,9  and  evaluating  Eq . 8 gives. 


E (z)  = /2a.  exp(-a  2)  . (10) 

IT  3d. 

This  equation  describes  the  near  field  directly  ad- 
jacent to  the  grating  waveguide.  Using  the  optical 
reciprocity,  theorem  we  apply  the  above  results,  ob- 
tained through  analyzing  the  decoupling  process,  to 
the  coupling  process  to  derive  the  coupling  efficiency 
(n^)  of  light  diffracted  into  a guided  mode. 

For  an  infinitely  long  output  grating  waveguide 
characterized  by  a , all  the  initially  guided  light  is 

3 

coupled  into  a radiation  mode  and  has  an  amplitude 
distribution  E^(z)  given  by  Eq.lO.  If  the  same  grating 
waveguide  is  used  as  an  input  coupler,  a 100%  coupling 
efficiency  can  be  obtained  if  the  input  amplitude  dis- 
tribution E\ (z)  has  the  field  distribution  E^(z)  and 
propagates  in  the  opposite  direction.  In  most  cases, 
E^(z)  i-  E^(z)  so  that  some  incident  energy  is  either 
reflected  or  transmitted.  Therefore  the  diffraction 
efficiency  for  light  coupled  into  the  waveguide  (n^) 
is  equal  to  that  fraction  of  the  incident  field  that 
corresponds  to  E^(z).  Therefore  , where  L is  the  coupler 
length , 


[/qEi  (z)E*  (z)dz]^ 

/n  E.E*dz/„  E E*dz 
Oil  0 r r 


(11) 


Because  the  decoupling  second  order  is  the  only  order 

illuminated  in  coupling  into  the  waveguide,  for  our  coupler 

the  factor  should  appear  in  Eg. 11.  If  a is 

d d d a 

constant  and  E^(z)  is  a uniform  plane  wave  with  length 
L then  Eq.ll  becomes, 

’’c  - dTTvT)  l5V>Il-exP(-o-aWl  .(12) 

'^d  ^ 


Eo..  12  indicates  has  a broad  maximum  with  a peak 
value  of  approximately  81%  for  a L = 1.25  and  zero 

3 


first  order  efficiency. 
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Calculation  of  Distributed  Diffraction  Efficiency 


By  TCP 

Alferness^  analyzed  the  wave  propagation  in  thick 

holographic  gratings  by  decomposing  a thick  grating 

into  thin  gratings.  Such  a method  can  be  extended  to 

analyze  the  optical  propagation  in  the  grating  inside 

the  waveguide  if  we  consider  the  mode  of  the  waveguide 

as  consisting  of  two  plane  waves  propagating  at  angle 

0 or  -0  [0  =sin  ^ (n  /n  )] . The  index  profile  of  the 
m mm  mg 

grating  inside  the  waveguide  can  be  expressed  as 

n(z)  = n (z) +n,  cos  (2iTf  (z-z  ) ) +n„cos  (4iTf  (z-z  ))  (13) 

g 1 X 2 X 

where  f is  the  frequency  of  the  fringes  in  any  constant 
X plane  (not  the  inverse  of  the  perpendicular 
distance  between  frinae  planes),  z^  = xtan(d))  where 
is  the  slant  angle  of  the  fringes,  and  n^  and  n^ 
are  the  fundamental  and  second  harmonic  index 
modulation,  respectively. 

We  examine  the  properties  of  the  optical  propa- 
cration  in  the  waveguide  of  thickness  D by  decomposing 
it  into  a series  of  thin  slabs  of  equal  thickness  Ax; 
where  the  magnitude  of  Ax  is  sufficiently  small  so 
that  each  slab  acts  as  a thin  grating.  The  amplitude 
transmittance  resulting  from  a plane  wave  incident  at 


81 


w 


0^  with  input  spatial  frequency 


n sin  (0  ) /\  is, 
g m 


T,^(z,x)  = exp  ( j 2nn  (z ) Ax/Xcos0  ) 

A m 

= exp{j2TTn  Ax/Xcos0^)  ? (j)*^J  (b^)  (14) 

9 q=-oo  ^ 

oo 

X exp  [ j27rqf  (z-z^)  ] Z (j)^J  (b2) 

X exp[  j4Trpf  (z-z^)  ]exp[j2n^Axcos0  /A] 
n g m 

where  >lq(b)  is  the  q order  Bessel  function, 

b.  = 27TAxn./Acos  0„  (for  i=l,2)  and  z = nAxtan  (i> ) . 
lira  n 

The  final  exponential  term  is  the  propagator  term  of  the 
incident  wave  over  the  slab  thickness  Ax. 

As  the  plane  wave  propagates  through  the  thickness 
Ax  in  the  waveguide,  it  diffractes  into  an  infinite 
set  of  plane  waves  of  spatial  frequency  (2p+q)f  and 
with  amplitudes  given  by  the  product  of  Bessel 
functions  of  order  q and  p.  Each  diffracted  v/ave 
then  becomes  the  incident  wave  of  the  next  thin  slab 
and  produces  another  set  of  diffracted  waves  according 
to  Eq.l4,  and  so  forth.  The  process  is  continued  until 
we  reach  the  other  boundary  of  the  waveguide,  where  the 
original  guided-wave  undergoes  total  internal  reflec- 
tion with  angle  -0  . Each  diffracted  wave  will  radiate 


J 


8^ 


into  the  air  or  substrate  region  at  an  angle  according 
to  Snell ' s law . 


Our  application  of  thin  grating  decomposition  to 
the  analysis  of  a waveguide  structure  containing  a 
slanted  grating  assumes  input  at  the  second  order 
Bragg  angle.  The  output  amplitudes  A(M)  resulting 
from  this  initial  condition  are  given  by 


A (M) 


(15) 


v^hich  is  evaluated  in  the  limit  as  Ax  -*■  0 and  m 
such  that  nAx  = 0,  and  where  m is  the  number  of  thin 
film  slices,  the  column  vector  is  the  initial  input  vector, 
and  H is  a 3x3  unitary  matrix  that  describes  the  coupling 
between  three  plane  waves  over  a thin  slab  of  thickness  Ax. 
We  are  then  required  to  perform  the  multiplication 
indicated  by  Eq.l5  with  initial  conditions  correspond- 
ing to  illumination  at  the  second  order  Bragg  angle. 

Referring  to  Fig. 3 we  define  C = cos0  and 

o m 

= cos0^.  By  inspection  of  Eq.l4  the  3X3  coupling 
matrix  for  each  thin  grating  becomes. 
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FIG.  3 

The  Bragg  Diagram  for  the  second  order  diffraction  in  a slanted  grating 
(represented  by  grating  vector  K)  which  is  constructed  by  interference 
of  pJane  waves  A and  D. 


J 


where , 


.^1 

^C~  ® 1 

1 

2ttAx  n 

>1  = [C^-CQ-Xf^tanc}.] 


01+  02 


f = ^(sin0, -sin0_) 
z A 1 2 

nn . Ax 

6.  = — (i  = 1,2; 

1 ^ 


with  0^  and  02  being  the  grating  construction  angles 
and  n^^  and  n2  being  the  fundamental  and  first  harmonic 
index  modulation.  To  facilitate  the  matrix  multiplica- 
tion of  Eq.l6  we  transform  the  coupling  matrix  into  a 

diagonal  matrix  H'  such  that  the  power  of  H becomes 

t h 

the  M power  of  the  diagonal  elements  of  H'.  The 
diagonal  elements  of  H'  (o^)  are  the  eigenvalues  of  H 
and  satisfy  the  matrix  equation, 


det(H-a^I)  = 0 . 
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For  convenience  we  adopt  the  form  of  the  eigenvalues 
as,  0^  = e ^*^i  (for  small  Ax)  . The  eigenvalue 

problem  is  to  find  solutions  which  satisfy 


^ . + p’lh  .-I 


' . +p'ii;.  + q'  ip . +r'  =0 
1 1 


(a  relationship  derived  from  Eq.l6  and  18) 


where , 


P'  = 

-^2  .2  .2 

q'  . ^ - !i_ 

V2  ■ 


The  solutions  of  Eq.l9  become. 


= U '+  V '+  P V3 


(U'+V*)  ^ (U-V)  P' 

~1~  ^ ~2~  ’^■3-  X 

-(u'+v*)  (u’-V)  /—  PJ. 

2 ~ 2 ~ 3 


where , 


v'.  rr-  ( 


with, 


a'  = 

i(3q’-p 

(26) 

b'  = 

1 /o  .3 

27  (2p 

'9p'q' 

+ 27r' ) . 

(27) 

The  output  amplitudes  of 

the  grating  in 

terms  of  the 

diagonal  matrix  b 

lasis  vectors 
ill) ' „ 

e-’^o  0 

become 

i 

0 

■ f 

■ M 

E(M)  = [H']’^E(0)  = 

0 

e-*  ^1 

0 

E(0) 

(28) 

0 

0 

e3^-'‘’^0 

0 

0 

= 

, 0 

0 

E(0) 

0 

0 

0 

0 

= 

0 

e=>'^0 

0 

E(0) 

(29) 

_ 0 

0 

e=>^^ 

The  associated  eigenvectors 

satisfying 

the  matrix 

equation. 

(H  - 

o .I)E 

. = 0 

(30) 

can  be  found  in  terms  of 

the  original  basis  vectors  to 

be, 


87 


where. 


(31) 


-(Y3>l^i-YiY2) 
Yj^Y3"Y2  (4'i+2C- 


C,  1/2 
^0 


(32) 


_ _C  (^.^2^.)  (r^^.-y^Y2)  C 

"i  - Cq  Yi[YiY3-Y2(V2^1^  S 


with. 


^1  = 


M6j^  TTn^D 


/CqCi  ^>^^0^1 


^2  = 


Yo  = 


= 


M62 

■nn^D 

/CqC2 

M(S^ 

irrij^D 

" 

/c  c 
•^1^2 

A/C^C2 

M4>, 

i 

2 

0 

i-> 

0 

0 

1 

(34) 


2TTDn 

c 
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ip^  in  Eq.32  and  33  are  obtained  from  Eq. 21-23 


using  Eq. 24-27  by  substituting  p,q,r  for  p',q',r' 


where , 


P = Mp'  = 25, 


w2  , 222 

q = M q'  = -Yj^  'Y2  -Y3 


r Mr  * 

Therefore  in  the  diagonal  representation  of  the 

I J 

coupling  matrix  H , the  total  amplitude  after  the  M 
thin  grating  is, 


E(M)  = 


e^^O 


where  the  input  column  vector  is  the  initial  condition 
in  terms  of  the  eigenvector  basis,  obtained  by  the 
transformation  relation  of  Eq.31  using  the  initial 
conditions  of  Eq.l5.  The  final  output  amplitudes  are 


obtained  by  transforming  the  output  E(M)  into  the  out- 
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put  A(M)  by  using  the  inverse  transformation  from 
Eq.30.  This  procedure  results  in  the  amplitude  of 
each  diffracted  order  as  being  given  by, 


^0  R, 


( 


(37) 


^2  = 


where , 


= 


^1-^0 


yo^i-UiVo 


^2‘^0 

- ( ■:■  - ^..  ..  ) 


(40) 


Rl  = Mq^v^-V2)  t Mi<V2-Vq)  + ^2 


(41) 


(42) 


The  amplitude  of  each  diffracted  order  as  given 
by  Eqs. 37-39  can  now  be  inserted  into  Eq.2  to  obtain 


the  diffraction  efficiency  of  each  order 
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r 


The  radiation  loss  coefficient  (a^)  can  then  be  deter- 
mined for  use  in  Eq.l2  giving  the  coupling  efficiency 
of  our  second  Bragg  order  waveguide  coupler. 
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Results  and  Discussion 


Using  the  approach  we  have  described  we  calcu- 
lated the  radiation  loss  coefficient  for  the  following 
four  representative  holographic  grating  couplers: 
a Slanted  or  unslanted  holographic  grating 
(with  same  periodicity  in  the  z-direction) 
in  the  waveguide  film  utilizing  the  first 
order  Bragg  effect. 

b Slanted  or  unslanted  holographic  grating  on 

top  of  the  waveguide  film  utilizing  the  first 
order  Bragg  effect  which  results  from  the 
interaction  between  the  evanescent  wave  of 
the  waveguide  and  the  holographic  grating, 
c Slanted  holographic  grating  in  the  waveguide 
film  utilizing  the  pure  double  diffraction 
process  (n2=0) , producing  the  second  order 
Bragg  effect. 

d Nonlinear  slanted  holographic  grating  in  the 
waveguide  film  utilizing  second  order  Bragg 
effect  (n2  0)  . 

In  case  b the  iiiteraction  product  of  the  evan- 
scent  wave  of  the  waveguide  and  the  holographic- 


grating  can  be  viewed  as  a plane  wave  incident  with 


mode  angle  0^  upon  the  holographic  grating  of  thickness 

Ax.  Since  the  guided  wave  is  incident  upon  a boundary 

of  the  waveguide,  the  reflected  wave  shifts  along  the 

propagation  direction  with  respect  to  the  position  that 

one  would  expect  from  geometrical  optics.  Such  a shift 

is  known  as  the  Goos-Hanchen  shift  (6z).  Thus  we 

should  express  the  depth  of  the  ray  penetration 

6z/ 

(evanscent  tail)  Ax  = 2tan6  . 

m 

The  diffraction  efficiencies  in  cases  a and  b 
were  calculated  with  the  numerical  analysis  (Eg. 14) 
and  cases  c and  d were  calculated  with  the  closed 
form  expressions  (Eq. 37-39).  Comparing  the  radiated 
loss  coefficient  of  the  four  generalized  coupler  con- 
figurations (Fig. 4-9)  it  is  apparent  that  the  slanted 
grating  of  case  a has  the  largest  radiated  loss  co- 
efficient (RLC)  with  the  slanted  coupler  of  case  d 
having  the  next  largest  RLC.  Approximately  the  same 
coefficient  results  for  the  unslanted  qrating  of  case 
a and  the  gratings  of  cases  b and  c . 

In  case  a (Fig. 4 and  5)  the  RLC  is  shown  to  be 
strongly  dependent  on  grating  slant  angle,  and 
the  coefficient  for  the  slanted  grating  is  two  orders 
of  magnitude  higher  than  that  of  the  unslanted  grating. 
We  note  that  in  thi s unslanted  grating  (Fig. 5)  that  the 
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FIG.  4 

Radiation  loss  coefficient  (um”^)  verses  n^  for  the  • 
D - 2.0pm  of  a case  a slanted  gratinq. 

.1. 
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TE^  radiation  loss  coefficient  (ym“^)  verses  n^^  with  n2  “ 0 for  the 
thickness  D = 2.0uni  of  a case  c slanted  grating. 


INDEX  MODULATION  (N,) 


FIG.  9 

TE^  radiation  loss  coefficient  verses  with 

for  the  thickness  D > 2.0um  of  a case  d slanted  grating 


TEj^  mode  has  a lower  loss  coefficient  than  either  the 
TEq  or  TE2  modes.  An  examination  of  the  electric 
field  pattern  of  the  TEj^  mode  reveals  a 180“  phase 
shift  and  therefore  a reduction  in  the  diffracted 
field  strength  as  the  wave  propagates  across  the  wave- 
guide. 

In  case  b the  radiation  loss  coefficient  is  ob- 
served to  have  a slight  dependence  on  the  grating 

slant  angle  (Fig. 6 and  1),  due  to  the  small  pene- 
tration depth  of  the  evanscent  tail  into  the  grating. 

If  the  grating  thickness  is  larger  than  the  evanescent 
tail  the  radiated  loss  coefficient  would  be  independent 
of  the  thickness. 

In  case  c,  (Fig. 8)  utilizing  the  double  dif- 
fraction process  the  second  order  Bragg  effect  coupler 
can  be  constructed  if  the  required  index  modulation 
can  be  achieved  inside  the  waveauide  film. 

In  case  d (Fia.9)  it  is  evident  that  an  effective 
second  order  Bragg  effect  coupler  could  be  constructed 

when  a nonlinearity  exists  in  the  film  (n2^  0) . 

4 

Kogelnik  and  Sosnowski  used  a thick  dichromated  gelatin 
holographic  phase  grating  to  achieve  70%  coupling  efficiency. 
Such  high  coupling  efficiency  can  be  attributed  to  the 
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first  order  Brsgg  diffraction  process  which  redis- 
tributes the  energy  from  various  diffracted  beeuns  into 
the  guided  beam.  Because  it  is  difficult  to  construct 
a waveguide  with  a sufficiently  high  spatial  frequency 
and  fringe  slant  to  provide  for  first  order  Bragg 
coupling  we  choose  to  employ  the  second  order  Bragg 
effect. 

Alferness^  has  demonstrated  that  the  second  order 
diffraction  process  can  result  from  both  double  dif- 
fraction through  the  fundamental  grating  frequency 
and  single  diffraction  throught  any  existing  second 
harmonic  grating  frequency.  If  the  second  harmonic 
component  is  negative  relative  to  the  fundamental, 
the  two  diffraction  processes  add,  giving  the  highest 
level  of  second  order  diffraction.  Utilizing  this 
second  order  Bragg  effect,  we  could  produce  in  an 
optical  waveguide  a Bragg  diffraction  grating  which 
couples  light  through  its  second  order  diffraction 
angle  into  a guided  mode  of  the  waveguide. 

A practical  and  efficient  grating  waveguide 
coupler  must;  i)  consist  of  a lossless  and  scatter-free 
material,  ii)  suppress  unwanted  diffraction  orders 
and  iii)  possess  a highly  spatially  modulated  re- 
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fractive  index  to  achieve  strong  waveguide  coupling, 
with  a thin  or  moderately  thick  (2-lOym)  layer  of 
phase  material  acting  as  a waveguide , our  grating  wave- 
guide could  meet  all  of  the  grating  coupler  constraints. 
This  coupler  would  be  fabricated  with  relative  ease 
because  of  the  low  grating  frequencies  required  for 
the  second  order  effects,  is  completely  integrated 
into  the  waveguide  material  and  has  the  potential  to 
achieve  high  coupling  efficiencies  from  air  to  wave- 
guide modes. 

In  summary,  a simple  geometrical  optics  concept 
is  presented  in  which  a grating  waveguide  coupler 
analysis  is  divided  into  independent  planar  waveguide 
and  diffraction  grating  problems.  The  analysis  of 
the  grating  is  accomplished  with  a closed  form  so- 
lution of  the  diffraction  efficiencies  derived  from 
the  thin  grating  decomposition  method.  Using  the 
conventional  analysis  of  a planar  waveguide  the 
distributed  radiation  loss  coefficient  of  any  grating 
waveguide  coupler  can  be  determined  from  the  separate 
grating  analyses.  The  results  of  our  analysis  are 
completely  compatible  with  those  of  the  analyses 

utilizing  the  perturbation  approach^  and  the  electro- 

2 

magnetic  boundary  value  approach  . However,  our 
analysis  has  the  advantages  of  being  conceptually 
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simple,  physically  intuitive,  with  results  obtained 
through  minimal  computation. 
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Abstract 

Algebraic  experssions  are  derived  for  the  diffraction 
efficiency  of  the  first,  second  and  third  orders  of  a holographic 
grating  recorded  in  a thick,  pure  phase  material  for  readout 
at  the  third-order  Bragg  angle.  The  results  demonstrate  the 
importance  of  triple  diffraction  by  the  fundamental  grating 
as  well  as  diffraction  by  harmonic  gratings.  The  method  can 
be  extended  to  analyze  the  first  order  Bragg  diffraction  at  the 
first  order  Bragg  angle  with  much  greater  accuracy.  It  is  shown 
that  strong  cross  coupling  can  be  obtained  from  triply- incoherently 
constructed  gratings. 
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In  recent  years,  the  thick  grating  has  been  studied  extensively 
and  a wide  range  of  applications  have  been  proposed,  such  as  highly 
efficient  diffraction  gratings,  narrow-band  optical  filters,  thlck- 
gratlng  optical  elements,  grating  deflectors  and  modulators,  grating 
couplers  In  Integrated  optics,  and  distributed-feedback  lasers  and 
distributed  Bragg-ref lector  lasers.  The  well-known  applications  of 
thick  holograms  Include  high  capacity  Information  storage,  color  ho- 
lography, and  whlte-llght-reconstructlon  holography.  In  essence,  the 
thick  hologram  can  be  regarded  as  recordings  of  a multiplicity  of 
thick  gratings. 

The  diffraction  of  a plane  wave  by  a thick  sinusoidal  grating  at 
or  near  Bragg  Incidence  was  studied  by  Kogelnlk  and  by  Burckhardt. 
Kogelnlk  obtained  a closed-form  expression  for  the  diffraction  efficiency 
at  the  first  order  Bragg  angle  by  employing  coupled-wave  theory.^ 
Burckhart  solved  the  exact  electromagnetic  boundary  value  problem  and 
obtained  numerical  results  by  using  the  computer  to  find  the  eigenvalues 
of  a matrix.  His  theory  considered  the  energy  exchange  among  higher 
orders  so  that  no  assumption  of  strong  higher  order  Bragg  suppression 
Is  required,  which  accoutns  for  It  being  used  In  gratings  of  arbitrary 
thickness.  However,  the  drawback  of  such  a method  lies  In  Its  dif- 
ficult Implementation  and  lengthy  computation.  Recently,  Alferness^ 
Introduced  the  thin  grating  decomposition  method  In  a single  grating 
case,  which  offers  conceptual  simplicity,  general  applicability  and 
ease  of  Implementation  In  the  computer.  His  method  has  been  used  In 
analyzing  second  order  diffraction  and  guiding  effects  In  thick  gratings. 
In  this  paper,  his  method  Is  extended  to  analysis  of  third  order 
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diffraction  of  thick,  high  index,  sinusoidally  modulated  gratings  which 
are  assumed  to  be  lossless.  The  surface  reflections  of  such  gratings 
can  be  neglected  in  the  derivation,  because,  in  practice,  they  can  be 


eliminated  by  anti-reflection  coatings.  Even  if  there  are  surface 
and  internal  reflections,  the  transmitted  diffraction  efficiency  can 
be  corrected  by  a multiplicative  transmittance  factor. 

Theory  of  Third  Order  Bragg  Diffraction 

Here  we  examine  a grating  constructed  interferometrically  in  a 
thick  emulsion,  as  shown  in  Fig.  1.  The  propagation  vectors  of  each 
of  the  recording  beams  are  in  the  xz  plane  and  are  symmetrical 
about  the  z-axis.  Therefore,  the  fringe  planes  are  perpendicular  to 
the  plate  and  there  is  no  variation  along  the  y axis.  The  thin 
grating  decomposition  method  is  valid  for  the  incident  readout  wave 
polarized  along  the  y axis.  When  we  consider  a pure  phase  material 
and  the  possibility  of  a film-generated  nonlinearity,  the  refractive 
index  inside  the  emulsion  can  be  written  as 

n(x)  = no  + nj  cos2irfx  + n2  cos4irfx  + no  cosbirfx  (1) 

where  ng  is  the  average  bulk  index  and  nj , n2  and  no  are  respectively 
the  index  modulations  of  the  fundamental,  second  and  third  harmonics 
of  spatial  frequency  f.  In  the  thin  grating  decomposition  method,  we 
assume  the  thin  vertical  slab  of  the  emulsion  to  be  of  thickness 
Az.  Since  the  fringe  planes  are  along  the  z direction,  the  thin 
slabs  are  identical.  The  spatial  phase  shift  caused  by  each  thin  slab 
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Figure  1.  Construction  geometry.  Pi  and  ?2  are  plane  vaves 

with  propagation  vectors  In  the  xz  plane  and  symmetrically 
Inclined  about  z axis.  The  fringe  planes  are  out  of  the 
paper  and  perpendicular  to  the  plane  of  the  plate.  The 
fundamental  grating  frequency  Is  f » 2(no/X)sln6  and  the 
emulsion  thickness  Is  D. 


incident  upon  by  a plane  wave  with  direction  cosine  cos0^  (angle 
measured  internally)  is: 

<l)(x)  = 2irn(x)  &z/Xcos6^  (2) 

where  X is  the  free-space  wavelength  of  the  readout  wave.  Since  the 
emulsion  is  a pure-phase  material,  the  amplitude  transmittance  of 
each  slab  can  be  written  as 


= exp(j(|)(x))  (3) 

The  amplitude  transmittance  of  each  thin-grating  slab  for  a wave 


incident  at  0^  can  be  expanded  as 

2irniAz 

(J>%tes“0~^  exp(j2.pfx)  X 
2irn2Az 

2im3Az 

(j)’'jr^XcoV07  ^ expCj&rrfx) 


where  we  have  dropped  an  unimportant  common  multiplicative  constant. 
The  total  field  emerging  from  the  first  thin  slab,  due  to  an  incident 

no 

plane  wave  of  spatial  frequency  f^  = sln0^  (the  component  of  spatial 
frequency  in  the  x direction),  is  the  product  of  the  input  amplitude 


A.  exp(j2irf  x)  and  the  amplitude  transmittance  of  the  slab: 
in  i 


out 


2TTni  Az 

^In  P=-“  ^^^'^'^p^Xcos©. 


+«? 


27rn2Az 


2trn3Az  2ff  ng  Xf 

^ (t S — ) exp  j— j — Az[l  - V-— 

r=-oo  r Xcos0  X ng 


expfj2Tr(f  + pf  + 2qf  + 3rf)x] 


(5) 
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where  we  include  the  propagator  of  the  incident  wave  over  the  slab  thick- 
ness Az. 

In  his  coupled-wave  analysis  for  readout  near  the  first  order 
Bragg  angle,  Kogelnlk  assumed  that,  for  gratings  with  a Q factor 
greater  than  10,  [Q  = ^-^f^D]  only  the  zeroth  (undiffracted)  and  one 
first  diffracted  order  are  important.  Bergsteln  and  Kermish^  have 
shown  that  in  a pure  sinusoidal  grating  with  4noni /(Xf )^<< 1,  only  the 
zeroth  and  mth  orders  are  significant  for  the  readout  at  the  mth  order 
Bragg  angle;  their  condition  implies  that,  in  addition  to  a large  Q 
factor,  a correspondingly  low  index  modulation  is  required  if  only 
two  orders  are  to  be  significant.  For  typical  gratings  or  special 
grating  profiles  (such  as  square,  triangular,  sawtooth  profile)  with 
a moderate  Q factor,  the  necessary  condition  for  only  two  orders 
does  not  apply  to  the  readout  at  the  higher  order  Bragg  angle. 

Therefore,  when  we  analyze  diffraction  at  the  third  Bragg  angle,  we 
have  to  include  four  significant  orders.  For  an  incident  wave  at 
z = 0 of  spatial  frequency  f q , the  spatial  frequency  of  the  first, 
second,  and  third  orders  are  fo  - f,  fo  - 2f,  and  fg  - 3f.  We 
write  the  amplitude  of  each  of  these  plane  waves  as  one  component 
of  a column  vector: 


Ac 

Al 

A2 

A3 


(6) 


where  the  subscripts  correspond  to  the  order  numbers.  The  result  of 


109 


each  thin-grating  slab  upon  the  propagating  waves  can  be  written  as  a 
coupling  matrix  which  describes  the  multiplicative  amplitude  and  phase 
change  effected  by  the  thin  grating  In  diffracting  the  jth  order  into 
the  ith.  The  1,  j matrix  elements  can  be  found  from  Eq.  5,  by 
choosing  p,  q,  and  r,  such  that, 

f^  = fj  + pf  + 2qf  + 3rf  (7) 

By  inspection  of  Eq.  5,  the  coupling  matrix  of  each  thin-grating 
slab  can  be  written  as 

, ^2j  ^34 

eJ«0  i— 

Cl  -^ci  -^cq 

6l-  ^l4  ^24 

gjai 

•Co  -’ci  ■’co 

624  ^l4  4 ^l4 

4_LJ“i  J“i  1— 

634  ^24  *^14  4 

1—^^“°  4__eJ“l  e^“° 

CO  -^co  -^ci 

where 


^1 

” irn^Az/X  1 =■  1,  2,  3 

(9) 

- 2rnoAz  cos6  /X  Z - 0,  1 

u 

(10) 

- cose^ 

(11) 

Since  we  can  make  Az  arbitrarily  small,  we  can  use  the  first  order 
approximations  of  the  Bessel  functions,  Jo(x)  = 1;  Ji(x)  = x/2; 
and  Jq(x)  = 0,  q > 2 for  x<<l,  and  Ji^ls  also  neglected.  We  have 
also  used  [1  - (—  - cos6«  and  J_i  (x)  » -J](x).  Since  the 

output  from  a thln-gratlng  slab  is  incident  upon  the  next  slab,  the 
total  amplitude  after  the  Nth  slab  Is 
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j 


\ = h\  (12) 

Because  the  readout  wave  Is  a single  plane  wave  of  unit  amplitude  and 
zero  phase,  we  have 


Here  we  analytically  carry  out  the  matrix  multiplication  of  Eq.  13 
of  the  readout  wave  Incident  at  the  third  Bragg  angle  In  the  limit 
N ^ Az  -*•  0,  such  that  NAz  = D,  the  emulsion  thickness.  The  matrix 
multiplication  of  Eq.  12  Is  greatly  simplified  by  transforming  It 
to  a basis  In  which  H Is  diagonalized.  The  advantage  of  a diagonal 
representation  Is  that  the  Nth  power  of  H becomes  simply  the  Nth  power 
of  each  of  the  diagonal  elements.  The  diagonal  elements  of  H in  this 
representation  are  the  eigenvalues  of  H and  the  basis  vectors  are  the 
associated  eigenvectors.  The  eigenvalues  of  H satslfy  the  matrix  equation 
det  (H  - 0^1)  - 0 (14) 

where  I Is  the  Identity  matrix.  It  Is  convenient  to  make  the  sub- 
stitution 

= -j  (1  - cJj^) 

The  solutions  of  Eq.  14  are: 

-A  ± /a2  + 4B 

’<'1.2  = 2 (15) 

-C  ± v^C^~r4D 

’"3.’* 2 (16) 


where 


6i  63 

A - (J)j 

Cl  CO 


- 6163  + 62  + ci634>j“  26j62 


(18) 


C 


4>1 


+ £1  + £1 

Cl  ^0 


(19) 


j)  , fil^  - <^163  + 62^  - Cl63(^^-^  26162 
cqCi 

2irno 

<I>1=  -J — (ci  - Cq) 


(20) 

(21) 


Since  we  can  make  Az  arbitrarily  small,  t|<j^«l,  therefore,  the  eigen- 
values are  = 1 - j\l>^  = exp(-jij)^).  The  associated  eigenvector 
satisfies  the  matrix  equation: 


(H  - o^DE^  = 0 

and  can  be  found  in  terms  of  the  original  basis. 


’Eo‘ 

[ 

El 

II 

E2 

_E3_ 

Ll 


(61  - 62)  ® 

(«i  - 62)  ® 

Cl  ('i'3  + 4£) 


(6]  - 62) 


Cl ('J'4 

‘Lole-J'f'l 


ci('f'l  - To) 

(61  - 62)  ® 

63 


ci(X-2- :ET) 


(61  - 62) 


(61  - 62) 
(61  - 62) 


(61  - 62) 


(23) 


From  the  Initial  condition  of  A (Eq.  13)  and  from  the  transformation 
of  Eq.  24,  the  initial  condition  in  the  E^  basis  is: 


E = 


(25) 
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Therefore,  in  the  diagonal  representation,  the  total  amplitude  after  the 
Nth  grating  Is 


N 

1 

re-wii 

p-J’J'2 

1 

®-jNi^2 

E(N)  - 

X 

1 

X 

1 

®-JN<|/4 

To  find  the  amplitude  of  the  first,  second  and  third  order  after  the 
Nth  thin  grating,  we  transform  back  to  the  original  basis  via  the 
inverse  of  the  transformation  in  Eq.  24,  thereby  obtaining: 


t(.3E3(N)  - KE2(N)  i(>2Eo(N)  " i(»lEi(N) 


(27) 


Ai(N) 


6l 

Cl 


i^l  - t(<2)(E2(N)  - E3(N))  - (i|^3  - i(;4)(Eo(N)  - EjCN)) 

2^  (i|/3  - - \l>2) 


■1 


(28) 


6i  :|^i^-(’i'l  - <^2)(E2(N)  - E3(N))  + (Jl>3  - i|^4)(Eo(N)  - Ej(N)) 
^2(N)  2^  (i/,3  - - \p2) 


A (N)  = -|[(- 


i{-3E3(N)  - KE2(N) 
(<l'3  - 'J'4) 


') 


i|j^0(N)  - 'l'lEi(N) 

^ ('J'2  - 'I'l) 


(30) 


with  NAz  = D,  we  define 


6jN  irniD 


yi  = 

, .1/2  = ,,  .1/2 
(cocj)  X(cqci) 

(31) 

62N 

Y2  “ 

(c.c.)'” 

(32) 

53N 

T3  • 

CO 

(33) 

Yl’  - 

Cl  CO 

(34) 

2imo 

^ - <(>iN  - — — D(ci  - cq)  (35) 


(29)  I 

J 

i 


113 


(36) 


r 


XI  - C - Yi'  + y3 

X2  * C + Yl'  ■ (37) 

Pi  * (xi  ^ + 4yii  P2  “ (X2^  + Ay  (38) 

Yll“  Vl  - Y2;  Y 22  “ Yl  + Y2  (39) 

XI '=  (xj  - 2ys);  X 2*  • (X2+  2y3)  (AO) 


Physically,  yi»  Y2»  Y3»  are  respectively,  the  total  grating  strength  of 

the  fundamental,  second  harmonic  and  third  harmonic,  and  C la  the  total 

phase  dispersion  (difference  In  propagation  factors)  between  the 

c. 

second  and  the  third  order.  The  diffraction  efficiency  • ~l^l  I ^ 
of  each  order  after  grating  thickness  D can  be  found  from  Eq.  27 
through  Eq.  30  by  using  Eqs.  15,  16,  and: 


'1,2 


In^os^l’ 

2r ' 2 2 


Yi^[( 


Pi 


sln^-^os^  sln^^os^  slnP|cos^ 

± -)  + ( -)  ± ( 7^ ^)]  (Al) 


P2 


Pi 


P2 


°3  ,0  “ i'  ^ sln^  cos^ln^  )] 

+ [ (^^In^sln^ --^^In^^ln^  ) + (cos^^os^  -cos^  cos^)  ] ^ (A2) 

We  readily  find  that 

DEo  + DEi  + DE2  + DE3  - 1 (43) 


Therefore,  the  energy  flow  In  the  z direction  Is  conserved — as  Indeed 
It  has  to  be,  since  for  a pure-phase  material  the  coupling  matrix  of 
Eq.  8 Is  unitary. 


PURE  SINUSOIDAL  GRATING  n2  - H3  - 0 


In  the  case  n2  “ n3  ■ 0,  all  the  energy  In  the  third  order  Bragg 


III+ 


diffraction  results  from  multiple  diffraction  through  the  fundamental 


grating.  An  interesting  property  of  the  third  order  Bragg  diffraction 
of  a pure  sinusoidal  grating  is  that  its  maximum  diffraction  efflency 
is  periodic  with  5.  By  inspection  of  Eq.  42,  we  can  express  the  con- 
ditions for  100%  DE3  as  either 

m • 1,  2,  3 . . . 

n = 1,  2,  3 . . . 

m i n (44) 

m = 1,  2,  3 . . . 

n = 1,  2,  3 . . • 

m i n (45) 

The  condition  imposed  on  K implies  that  for  a fixed  thickness  and 
readout  wavelength,  the  third  order  diffraction  efficiency  can  be 
100%  only  for  discrete  grating  frequencies.  In  Fig.  2,  we  plot  DE3  vs 
Yi  for  several  values  of  (cq/ci)^^^,  with  C “ The  maximum  value 

of  DE3  can  be  made  100%  for  proper  choice  of  Yi  and  the  ratio  (cq/cj)^^^. 


Pl/2  ■ nm 
P2/2  “ (m  + 1)tt 
rj '/2  = n 

C • (2m  + l)ir/2n 
or 

p j / 2 “ mu 
P2/2  = (m  + 2)it 
ri’/2  - (2n  + 1) 

C = 4(m  + l)Tf/(2n  + 1) 


1 


j 


i 
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Figure  2.  Diffraction  efficiency  of  the  third  order  at  the 
third-order  Bragg  angle  due  to  triple  diffraction 
(n2  “ 03  » 0)  versus  the  fundamental  grating  strength 
Yi  ” irniD/X  f co/q)  1/2  ^Ith  C * 5tr/2  for  several  values 
of  (co/ci)l/2. 


GENERAL  CASE:  n 7.  n ^ 0 


In  most  cases  of  interferometrically  constructed  gratings  of 
pure  phase  materials,  the  film  response  Is  nonlinear;  therefore 
we  must  Include  all  the  terms  nj,  n2»  and  n3-  The  necessary 
conditions  for  100%  diffraction  efficiency  In  the  third  order 
are : 


xi”  = 2nir  = Yi’  + Y3 


pj/2  = mu;  P2/2  = (m  + 1)it 


(2m  + - 4Y1Y2 

[2nu  - 2Y3] 


m > n 


(46) 


Because  the  condition  on  5 Include  Yi,  Y2  andY3»  DE3  can  be  as 
high  as  100%  for  a fixed  thickness  and  wavelength  and  for  any 
arbitrary  choice  of  spatial  frequencies  by  proper  choice  of  nj, 
n2  and  n3.  In  the  case  where  n^,  n2  are  each  zero,  the  result  Is 
Identical  to  the  coupled  wave  result  at  the  first  Bragg  Angle. 
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RESULTS  FOR  Yi/C.  V VC  « 1 

The  general  results  of  Eqs.  40  and  42  can  be  simplified  under 
the  assumption  of  index  modulation  small  relative  to  (Xf)^,  or  yi/C» 
Y2/C  <<!•  The  following  small  parameters  are  defined  as: 


ti  “Yi/C  = 2nonj/ (Xf ) 2 <<  l (47.) 

tz'Yz^?  = 2noni/(Xf)2  « 1 (48) 

^3  ='Y3/5  = 2noni/(Xf)2  <<  1 (49) 

ti'=  Yl'/C  = /^  2noni/(Xf)  2 « 1 (50) 

We  see  that  DE j from  Eq.  42  is  simplified  to: 

DE3  = sin2[Y3  - (51) 


The  three  terms  in  the  bracket  can  be  identified  as  three  possible 
modes  for  achieving  third  order  diffraction  through  the  three  Fourier 
components  f,  2f,  3f;  they  are  (1)  triple  diffraction  through  first 
harmonics,  (2)  combinations  of  diffraction  through  second  harmonics 
and  first  harmonics,  (3)  diffraction  through  third  harmonics. 

Eq,  51  shows  the  opposing  effects  of  mode  2 against  mode  1 and  3, 
if  n2  and  n3 are  positive  relative  to  ni.  Such  results  agree 
with  those  of  D.  L.  Jaggard®  in  his  extended  coupled  wave  analysis  of 
higher  order  Bragg  coupling  in  a periodic  medium;  also  with  those  of 
Su  and  Gaylord^  on  higher  order  diffraction  hyutllizing  the  Chu  and 
Tamlr®  rigorous  modal  approach. 
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EXTENSION  OF  ANALYTIC  RESULTS  FOR  THE  FIRST-ORDER  READOUT; 

We  can  increase  the  accuracy  of  the  analytical  results  for 
readout  at  the  minus  first  Bragg  Angle  by  including  also  the 
positive  first  and  minus  second  orders.  This  is  especially 
true  for  gratings  of  moderate  Q factor  or  very  high  index  modulation. 

At  the  first  Bragg  angle,  the  dephaslng  term  associated  with  the 
conjugate  first  order  is  the  same  as  that  for  the  second  order,  and 
in  general,  these  orders  can  be  of  comparable  Importance.  The 
coupling  matrix  for  the  first  order  readout  is  the  same  as  that  for  the 
third  order  readout,  except  for  the  modification  of  the  initial  condition: 

(52) 

The  minus  first  diffraction  order  of  the  results  from  the  general  four- 
waves  analysis  can  be  calculated.  From  the  results  we  can  see  that 
as  the  Q factor  Increases,  the  first  order  diffraction  efficiency 
gets  closer  to  100%  when  Yj  * if/2,  which  means  that  as  the  Q factor 
Increases,  the  higher  orders  can  be  neglected.  The  closed  form 
solutions  for  diffraction  efficiency  of  first  positive,  zeroth,  minus 
first,  and  minus  second  orders  can  be  obtained  by  changing  the  initial 
condition  to  those  of  Eq.  52. 

The  diffraction  efficiency  of  such  diffraction  orders  are: 
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Pi 


P 2 2 ^ Pi  ^ ^ 


(53) 


DEo,-i“  i ^ 

+ (cos^cos^  ± COS^-^OS^  )]2 

+ [(^^In^^os^  ± -psln^os^  ) 

Pi  2 2 P2  ^ 2 

+ (cos^sin^  - cos^in^  )]2  ^ 

With  small  C > a grating  can  be  considered  as  a thin  grating: 
there  is  no  Bragg  effect  on  the  diffraction  efficiency  and  the  maximum 
value  of  DEi  = 33%  agrees  very  well  with  the  first  maximum  of  Jj^. 


With  an  intermediate  C»  the  diffraction  efficiency  is  strongly 
dependent  on  C and  a rapid  monotonic  increase  in  DE_j  as  the  other 
orders  are  suppressed.  In  cases  of  very  large  DE_j  does  not 
depend  entirely  on  higher  harmonic  grating  frequencies.  The  dependence 
of  diffraction  efficiency  on  higher  harmonic  grating  frequencies  is 
quite  significant  as  we  can  see  from  Figure  3,  where  Yi  = 1.0, 

Y2  ” -0.5, and  various  Y3.  For  large  the  dependence  of 

the  diffraction  efficiency  on  third  harmonic  grating  frequencies  is 

reduced  drastically — almost  Independent.  From  Figure  4,  where  Yi  = 1.0, 

Y3  “ 0,  and  various  Y2»  see  that  even  for  large  5,  at  least  a 10%  increase 
in  diffraction  efficiency  can  be  expected  when  the  value  of  Y2  Increases 
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Figure  3.  First  order  diffraction  efficiency  at  the  first 
Bragg  angle  versus  Y3  with  a fixed  Yi  “ 1*  and 
Y 2 ’ for  several  values  of  C* 


Figure  4.  First  order  diffraction  efflcinecy  at  the  first 
Bragg  angle  versus  y 2 ® fixed  yi  “ !• 

Y3  =•  0.  for  several  values  of  C» 


16 
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from  0 to  0.5.  Therefore,  it  is  evident  that  even  for  the  first  order 
diffraction  it  is  necessary  to  include  the  higher  order  grating  com- 
ponents in  order  to  obtain  accurate  results. 

ANALOGY  WITH  MULTIPLE  GRATING  g - 0 

Figure  5 is  a construction  geometry  for  triply  Incoherently 
superimposed  gratings.  We  first  interfere  plane  waves  Pj  and  P 2, 
then  ?2  and  P3,  and  finally  P3  and  P4.  The  incident  wave  at  the 
Bragg  angle  of  the  first  grating  is  strongly  diffracted  by  it. 

This  diffracted  wave  is  at  the  Bragg  angle  of  the  second  grating  and 
is  strongly  diffracted  again.  And  it  finally  is  strongly  diffracted 
by  the  third  grating.  The  cross-coupled  order  via  such  triple 
diffractions  is  the  same  as  the  third  order  diffraction  in  the  single 
grating  of  "2  ” "3  ” 0*  proved  that  the  coupling  matrix  for 

the  multiple  grating  described  above  is  identical  to  Eq.  8,  except 
that  C “ 0.  The  result  for  DE3  will  be  the  triple  diffraction  for 
a triple  Incoherently  superimposed  grating  as  Illustrated  in  Fig.  5. 
The  triple  diffraction  results  are  shown  in  Fig.  6 for  different 
ratio  of 
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Figure  5.  Construction  geometry  for  multiple  gratings.  For  the 
incoherent  superposition,  plane  waves  Pj  and  ?2  are 
interfered  first,  then  P2  and  P3,  and  finally  P3  and 
P4. 
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